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FOREWORD 


This  report  was  prepared  by  the  McDonnell  Douglas  Corporation, 
McDonnell  Aircraft  Company,  Avionics  Systems  Technology  Department  under 
USAF  Contract  F33615-70-C-1333*  The  program  was  administered  by  the  Air 
Force  Avionics  Laboratory  (AFSC )  Systems  Avionics  Division,  Wright-Patterson 
Air  Force  Base,  Ohio.  The  Air  Force  project  engineer  directing  the  techni¬ 
cal  aspects  of  the  study  was  Herbert  Bartman,  AFAL/AAA. 

This  document  is  Technical  Report  AFAL-TR-71-155,  Part  I,  and  covers 
work  accomplished  from  January  1970  to  April  1971.  Part  II  will  describe 
verification  methods  and  results,  and  Part  III  will  be  a  users'  manual. 

This  technical  report  has  been  reviewed  and  is  approved. 


WILLIAM  A.  STUDABAKER 
Lt  Colonel,  USAF 
Chief 

System  Avionics  Division 


ABSTRACT 


A  computer  program  was  developed  to  predict  and  analyze  electromagnetic- 
interference  between  avionics  systems  on  aerospace  vehicles.  It  is  an 
analytic  tool  to  aid  the  user  in  establishing  intra -vehicle  electromagnetic 
compatibility.  The  program  is  used  to  derive  information  which  is  other¬ 
wise  impractical  or  difficult  to  obtain.  This  part  covers  the  investi¬ 
gation,  research,  and  development  phases  of  this  program.  A  study  of 
physical  principles,  previous  work,  and  practical  experience  identified 
four  predominant  interference  coupling  paths: 

0  Antenna -to-antenna 

0  Wire -to -wire 

°  Electromagnetic  field  (antenna) -to -wire 
°  Box -to -box 

Basic  equations  and  analysis  techniques  were  formulated  for  each  coupling 
path.  These  were  then  assembled  into  four  practical  computer  analysis 
programs  using  FORTRAN  IV  language.  The  programs  are  modular  and  thus  can 
be  run  separately  or  combined  into  one  integral  program.  This  report 
contains  a  complete  description  of  each  program  including  the  mathematical 
models,  basic  equations,  and  analysis  procedures  used.  Examples  of  the 
application  of  the  programs  together  with  sample  computer  outputs  are 
presented.  The  programs  were  developed  for  the  CDC  6600  computer,  but  they 
are  readily  ad.  ptable  to  other  machines. 
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SECTION  1 


INTRODUCTION  AND  SUMMARY 


Modem  aerospace  vehicles  contain  increasing  numbers  of  complex  and 
sophisticated  avionics  systems  that  must  operate  harmoniously  to  achieve 
mission  objectives.  Achieving  electromagnetic  compatibility  (EMC)  for  such 
avionics  is  a  complex  and  time  consuming  problem.  To  assist  the  engineer 
(whether  DCD  agency,  prime  contractor  or  sub-contractor)  in  establishing 
and  maintaining  intra -vehicle  electromagnetic  compatibility  (IVEMC',  the 
power  and  speed  of  modem  computers  can  be  utilized  with  si  gnificaut  1  ine 
and  cost  savings.  The  engineer  retains  responsibility  for  the  appropriate 
decisions  as  to  what  interference  specifications  and  control  measures  are 
required, but  with  the  tedious  computations  performed  by  the  computing 
equipment,  he  can  devote  more  of  his  time  to  the  design  decisions. 

McDonnell  Aircraft  Company  (MCAIR)  under  contract  to  the  Air  Force 
Avionics  Laboratory*  has  developed  analytical  techniques  for  predicting 
electromagnetic  interference  (EMI)  between  avionic  systems  on  aerospace 
vehicles.  The  techniques  derived  are  specifically  for  analyzing  coupling 
and  interference  propagation  paths.  The  result  is  a  series  of  interrelated 
computer  programs  to  perform  the  operations  indicated  by  the  mathematical 
models  developed  for  the  solution  of  IVEMC  situations.  The  development  of 
the  computer  programs  is  divided  into  three  phases:  investigation, 
development,  and  verification.  This  part  addresses  the  first  two  phases. 

The  third  phase  will  be  described  in  part  II. 

Investigation  Phase  -  During  the  investigation  phase,  the  basic  con¬ 
cepts  and  solution  approaches  were  established.  The  IVEMC  problem  was 
examined,  and  the  major  EMI  coupling  paths  were  identified  by  reviewing 
the  physical  principles  involved,  previous  investigation  results,  and 
McDonnell  Douglas  Corporation  (MDC)  and  USAF  experience  with  IVEMC  problems. 
The  basic  equations  and  techniques  for  the  mathematical  models  were  then 
established  and  verified  for  a  number  of  cases. 

Development  Phase  -  During  the  development  phase,  the  basic  equations 
and  techniques  were  assembled  into  practical  computer  programs.  The  input 
data  and  the  required  outputs  were  determined,  and  the  algorithms  were 
formulated.  The  algorithms  were  coded  for  the  CDC  6600  computer,  and  the 
programs  debugged. 

Verification  Fhase  -  A  number  of  tests  were  performed  to  ensure  the 
validity  and  accuracy  of  the  prediction  models.  These  included  screenroom 
component  tests,  full  scale  mock-up  tests  in  an  anechoic  chamber,  and  actual 
aircraft  tests.  Results  from  these  and  current  tests  will  be  included  in 
part  II  which  will  be  submitted  in  September  1971* 
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SECTION  2 


APPROACH  TO  PROBLEM 


2.1  ANALYTICAL  APPROACH 

To  design  for  intra-vehicle  electromagnetic  compatibility,  the  engineer 
must  have  the  ability  to  accurately  predict  levels  of  expected  interference. 
The  computer  programs  described  herein  are  designed  to  assist  him  in  making 
these  predictions.  The  engineer,  however,  must  assess  application  validity, 
since  these  programs  are  merely  models  of  the  actual  physical  situation. 

Of  necessity,  they  contain  many  assumptions  and  approximations  which  the 
engineer  must  take  into  account. 

These  programs  have  8  broad  capability  to  analyze  the  actual  case  at 
hand;  they  can  largely  reduce  reliance  solely  on  past  experience  and 
"rules -of -thumb."  The  speed  and  accuracy  of  the  computer  can  be  utilized 
to  obtain  a  picture  of  what  interference  will  occur,  how  severe  it  will  be, 
what  1  he  predominant  coupling  mode  is,  what  the  most  offensive  sources  are, 
ar. i  what  the  most  susceptible  receptors  are.  Knowing  this,  the  engineer  can 
make  decisions  as  to  what  can  be  done  to  reduce  or  eliminate  the  interference. 
The  computer  can  then  be  used  to  verify  the  adequacy  of  these  decisions. 

The  adequacy  of  applicable  EMI  specifications  can  also  be  determined,  and 
the  specifications  can  be  modified  as  needed. 

These  programs  meet  the  needs  in  establishing  IVEMC  on  multi-purpose 
air  vehicles.  Equipment  on  these  vehicles  must  be  interchangeable  from  one 
system  to  another,  and  often  identical  line -replaceable -units  are  used  on 
widely  different  aircraft  types.  Further,  these  vehicles  are  subject  to 
extensive  modification,  often  to  applications  far  different  from  that  of  the 
initial  design.  To  apply  over  the  broad  range  of  air  vehicle  EMI  analysis 
requirements,  these  programs  are  broad  in  scope  and  flexible.  They  will 
handle  the  large  number  of  EMI  situations  encountered  on  these  vehicles,  and 
they  can  he  easily  modified  for  additional  applications.  The  programs  can 
be  run  as  a  group  to  provide  total  system  SMI  data,  or  they  can  be  run 
Independently  for  a  particular  coupling  mode.  Thus,  if  antenna -to -antenna 
coupling  is  of  interest,  only  that  program  need  be  run. 

These  programs  can  also  be  used  to  develop  special  purpose  specifica¬ 
tions  for  limited  production  vehicles  such  as  spacecraft.  For  a  defined 
system  configuration,  the  subsystem  parameters  can  be  fed  into  the  computer, 
and  the  interference  levels  calculated.  Specifications  can  then  be  written, 
or  existing  specifications  can  be  modified  as  necessary  for  IVEMC. 

2.2  PROBLEM  DEFINITION 

The  first  step  in  establishing  these  IVEMC  analysis  programs  was  to 
determine  the  predominant  EMI  coupling  modes  and  paths.  Once  these  were 
determined,  mathematical  models  were  devised  for  each  mode.  MCAIR  had 
previously  developed  an  antenna -to-antenna  program  (Ref  1) ,  and  a  survey 
was  initiated  to  determine  other  coupling  modes  for  which  useful  prediction 
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programs  could  be  devie  d  within  the  scope  of  this  contract.  The  range  of 
USAF  interference  problems  was  examined.  Also,  an  extensive  literature 
search  was  conducted  to  review  past  approaches  and  limits  of  current  practice. 
Future  trends  in  avionics  were  surveyed.  From  these  surveys,  four  coupling 
paths  were  determined  most  predominant: 

°  Antenna -to -antenna  (transraitter-to-receiver)  coupling 

°  Wire -to -wire  coupling 

°  Electromagnetic  field -to-wire  coupling 

°  Box-to-box  magnetic  coupling 

Antenna-to-Antenna  Coupling  -  One  of  the  most  common  interference 
problems  on  aerospace  vehicles  results  from  coupling  from  a  transmitter 
antenna  to  a  receiver  antenna.  Large  numbers  of  radiating  and  receiving 
avionics  are  being  carried  on  modern  aircraft  with  a  trend  to  higher 
powered  transmitters  and  more  sensitive  receivers.  One  version  of  the  F-4, 
for  example,  has  transmitting/ receiving  subsystems  which  use  more  than  30 
antennas.  When  a  new  subsystem,  such  as  an  SCM  pod,  is  added  to  such  a 
vehicle,  determining  the  probable  interference  without  the  aid  of  a  computer 
is  a  large,  formidable  task. 

Wire -to-Wire  Coupling  -  Coupling  from  one  wire  to  another  within  a 
bundle  has  long  been  recognized  as  a  serious  source  of  SMI.  With  current 
trends  toward  higher  density  packaging,  more  compact  bundles  with  smaller 
wires  and  thinner  insulation,  the  problem  is  becoming  more  acute.  MCAIR 
performed  comprehensive  investigations  in  this  area  and  conducted  an 
extensive  literature  search  so  as  not  to  duplicate  previous  efforts. 

Electromagnetic  Field -to-Wire  Coupling  -  Experience  with  avionics 
systems  in  aircraft,  missiles,  and  spacecraft  has  shown  coupling  to  wire 
bundles  within  the  vehicle  from  external  RF  energy  to  be  a  major  EMC  concern. 
Usually,  the  energy  was  from  an  on-board  radiating  antenna  and  was  coupled 
onto  the  wires  via  gaps  or  openings  in  the  vehicle  skin.  Such  apertures 
include  radomes,  canopies,  landing  gear  doors,  access  doors,  and  engine 
intakes. 

Box-to-Box  Magnetic  Coupling  -  Another  common  problem  occurs  when  the 
lOO  Hz  magnetic  field  of  an  avionics  box  is  coupled  into  sensitive  circuits 
in  other  boxes.  The  usual  source  of  these  fields  is  from  power  transformers, 
but  other  high  current  devices  such  as  solenoids,  relays,  and  motors  can 
also  produce  them.  Common  receptors  for  such  interference  include  low  level 
audio  and  video  circuits  as  well  as  digital  logic  and  servo  circuits. 

Other  combinations  of  paths  and  coupling  modes  can  and  do  cause 
ird  '-i'  -v  no  ,  but  th  four  li  t  d  abov  arc  tb  •  bn.  ic  on'’  ;  .nd 
they  meet  the  criterion  of  providing  the  engineer  with  practical,  useful 
tools  for  establishing  IVEMC.  Further,  since  they  are  basic,  the  resulting 
mathematical  models  lend  themselves  to  widespread  applications  in  analyzing 
other  coupling  modes. 
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,jith  the  exception  of  the  wire-to-wire  program,  only  first-order 
effects  were  considered.  (See  Section  U.2. 1.)  In  the  past,  many  attempts 
have  been  made  to  analyze  multi -order  effects,  such  as  intermodulation. 

While  this  is  a  potential  source  of  EMI,  experience  has  shown  it  to  be  rare 
when  compared  to  direct,  first-order  interference.  Multi-order  EMI  usually 
occurs  only  in  special  situations,  as  for  example  two  transmitters  tuned  to 
a  specific  frequency  differential  simultaneously.  This  project  was  there¬ 
fore  directed  to  modeling  only  first-order  effects. 

P.3  UNDERLYING  PRINCIPLES 

For  maximum  usefulness,  the  analytic  tools  must  be  oriented  toward  the 
user  -  not  the  computer.  That  is,  the  computer  should  work  for  the  user, 
not  vice  versa. 

2. 3. 1  Data  Format 

These  programs  require  and  accept  data  which  is  normally  available  in 
several  forms.  The  format  and  codes  are  easily  learned  so  that  extensive 
training  and  computer  knowledge  is  not  required  of  the  user.  Also,  the 
data  formats  are  common  between  the  programs.  The  same  data  cards  describing 
a  wire  bundle,  for  example,  are  accepted  by  both  the  field -to-wire  and  wire - 
to -wire  analysis  programs.  This  way,  only  one  set  of  cards  is  needed  for  a 
given  set  of  data.  The  cards  are  merely  assembled  in  different  ways  for  the 
different  programs. 

?.  3 .  2  Default  Features 

Since  these  programs  are  intended  for  use  early  in  the  conceptual  and 
design  phases  of  vehicle  development,  they  have  to  be  usable  before  many  of 
the  basic  equipment  parameters  are  known.  In  this  event,  the  programs  have 
built-in  default  assumptions  based  on  the  applicable  military  specification 
or  on  mathematical  expressions.  Thus,  an  analysis  can  be  performed  with  the 
default  values  so  that  the  major  EMI  problem  areas  can  be  determined. 

Later,  when  the  actual  data  and  specifications  become  available,  new 
analyses  can  be  made  to  update  the  previous  ones. 

2.3.3  Program  Outputs 

After  the  data  is  read,  each  program  prints  a  summary  of  this  data 
including  any  default  values  that  were  inserted.  Thus,  a  record  of  the  data 
on  which  the  analysis  was  based  is  provided.  The  default  values  are 
identified  in  the  output  summary. 

The  analysis  results  are  provided  as  EMI  margins  in  decibels  (dB). 

(The  EMI  margin  is  defined  as  the  ratio  of  the  received  signal  at  a  given 
receptor  from  a  given  interference  source  to  the  susceptibility  level  of 
that  receptor.)  Thus,  the  degree  of  interference  can  be  readily  ascertained. 
Information  on  the  cause  of  the  interference  and  the  analysis  method  utilized 
are  also  given  in  the  output  for  each  interference  situation.  In  the  antenna 
to -antenna  program,  for  example,  the  maximum  EMI  frequency  and  the  maximum 
and  minimum  frequencies  of  coincidence  between  the  transmitter  and  receiver 
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are  given.  Where  applicable,  outputs  are  provided  that  are  directly 
relatable  to  military  specifications.  For  example,  the  wire-to-wire  pro¬ 
gram  gives  the  received  signal  and  susceptibility  levels  in  dBg  V/MHz  as 
used  in  the  EMC  specifications. 

The  programs  also  provide  extensive  diagnostic  and  error  data.  For 
example,  if  an  input  code  is  not  recognized  by  the  computer,  it  prints  out 
the  type  of  code  as  well  as  the  bad  data.  This  aids  the  user  In  locating 
a  mispunched  or  cut  of  sequence  data  card. 

The  programs  provide  for  the  use  of  classified  data.  The  classification 
is  inputed  for  each  subsystem,  and  the  highest  classification  is  printed  on 
the  top  and  bottom  of  each  page  of  output. 

2.3.1*  Mathematical  Model  Formulation 

The  mathematical  models  were  formulated  to  rapidly  survey  a  large 
amount  of  data  and  give  a  reasonable  indication  of  the  interference  that 
will  occur.  The  models  account  for  major  effects  on  interference  coupling 
(e.g. ,  aircraft  fuselage  shading).  The  analysis  results,  however,  are 
sufficiently  accurate  to  provide  for  correct  EMI  engineering  decisions. 

Thus,  a  balance  is  obtained  between  speed  and  accuracy. 

The  equations  and  computational  methods  that  comprise  the  math  model 
were  formulated  from  basic  electromagnetic  relationships.  These  relation¬ 
ships  Include  the  basic  definitions  of  capacitance  and  inductance, 

Maxwell's  equations,  and  the  wave  equation.  Where  approximations  were 
used,  care  was  taken  that  the  accuracy  was  retained.  Further,  the  models 
were  designed  to  give  answers  consistent  with  past  experience  and  "real 
world"  situations. 

To  ensure  that  all  potential  EMI  is  detected,  worst -case  assumptions 
are  used  for  vague  or  undefinable  parameters,  such  as  wire  separation 
within  a  bundle.  In  many  cases  statistical  limits  could  be  ’used,  but  in 
a  vehicle  development  phase,  where  these  programs  are  of  most  value, 
worst -case  analysis  is  desired. 
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are  given.  Where  applicable,  outputs  are  provided  that  are  directly 
relatable  to  military  specifications.  For  example,  the  wire -to -wire  pro¬ 
gram  gives  the  received  signal  and  susceptibility  levels  in  dB/i  V/MHz  as 
used  in  the  EMC  specifications. 

The  programs  also  provide  extensive  diagnostic  and  error  data.  For 
example,  if  an  input  code  is  not  recognized  by  the  computer,  it  prints  out 
the  type  of  code  as  veil  as  the  bad  data.  This  aids  the  user  in  locating 
a  mlspunched  or  out  of  sequence  data  card. 

The  programs  provide  for  the  use  of  classified  data.  The  classification 
is  inputed  for  each  subsystem,  and  the  highest  classification  is  printed  on 
the  top  and  bottom  of  each  page  of  output. 

2.3.4  Mathematical  Model  Formulation 

The  mathematical  models  were  formulated  to  rapidly  survey  a  large 
amount  of  data  and  give  a  reasonable  indication  of  the  interference  that 
will  occur.  The  models  account  for  major  effects  on  Interference  coupling 
(e.g. ,  aircraft  fuselage  shading).  The  analysis  results,  however,  are 
sufficiently  accurate  to  provide  for  correct  EMI  engineering  decisions. 

Thus,  a  balance  is  obtained  between  speed  and  accuracy. 

The  equations  and  computational  methods  that  comprise  the  math  model 
were  formulated  from  basic  electromagnetic  relationships.  These  relation¬ 
ships  Include  the  basic  definitions  of  capacitance  and  Inductance, 

Maxwell's  equations,  and  the  wave  equation.  Where  approximations  were 
used,  care  was  taken  that  the  accuracy  was  retained.  Further,  the  models 
were  designed  to  give  answers  consistent  with  past  experience  and  "real 
world"  situations. 

To  ensure  that  all  potential  EMI  is  detected,  worst -case  assumptions 
are  used  for  vague  or  undefinable  parameters,  such  as  wire  separation 
within  a  bundle.  In  many  cases  statistical  limits  could  be  used,  but  in 
a  vehicle  development  phase,  where  these  programs  are  of  most  value, 
worst -case  analysis  is  desired. 
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The  EM  is  computed  by  expanding  Equation  (3*1-)  into 


EM  P  +  TFS+SF  CLx-CLr-SR 


where 


P 

TFS 

SF 

CL 

X 

CL 


R 


the  transmitter  power  in  dBm. 
free  space  transmission  factor  in  dB. 
shading  factor  in  dB. 
transmitter  cable  loss, 
receiver  cable  loss. 


(3.2) 


The  term  "shading"  refers  to  signal  attenuation  due  to  propagation  around 
the  aircraft  exterior  surfaces  (fuselage  and  wings).  Computation  of  the 
terms  of  Equation  (3*2)  is  as  follows. 


Receiver  Model 


3R,  the  receiver  threshold  function,  is  represented  by  the  curve  shown 
in  Figure  1.  In  the  figure,  Ftr  and  Fcr  are  the  lower  and  upper 
frequencies  of  the  receiver  tuning  range,  and  Sm  is  the  minimum  sensitivity 
level  in  dBm.  The  selectivity  skirt  levels  are  represented  by  3Lp  (in  dB) 
”nd  skirt  frequencies  Far  and  Fjy,  which  are  obtained  with  the  receiver 
tuned  to  each  end  of  its  operating  frequency  range. 


Relative 
Power  Level 

(dBm) 


Log10(F) 
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FIGURE  1  RECEIVER  THRESHOLD  LEVEL  vs  FREQUENCY 


3  Transmitter  Model 


In  a  similar  manner,  the  transmitter  power  spectrum,  P,  is  represented 
as  shown,  in  Figure  2.  Only  parameters  for  the  fundamental  curve  shape 
are  specified  since  those  of  the  harmonics  are  assumed  identical  except  for 
downward  displacement  attenuation.  The  limits  of  the  operating  frequency 
range  are  specified  by  Fbtf  and  Fctf.  Skirts  are  represented  by  the  skirt 
L  *vel,  3L<r,  and  frequencies  Fa^f  and  F^tf*  ?m  is  the  peak  output  power  in 
dBm.  The  harmonic  curves  are  displaced  from  the  fundamental  by  a  factor, 
HF.  as  shown. 
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FIGURE  2  TRANSMITTER  EMISSION  POWER  SPECTRUM 


3.?.^  Propagation  Model 

The  free  space  transmission  factor,  TFS,  in  Equation  (3«?)  in  found  Ly 
using  the  Friis  transmission  equation  (Ref  2) 


TFS  Gx  +  Gr  +  20  log10 


X 

4  7T  D 


where 

Gx  gain  of  the  transmitter  antenna  in  dB 
Gr  gain  of  the  receiver  antenna  in  dB 
D  distance  between  the  two  antennas  in  meters 
X  the  wavelength  in  meters 


3 .  t .  1  Antenna  Gains 

Gy  and  Gp  are  generally  functions  of  direction  and  ar>j  determined  by 
ieflning  the  angles  0  and^,  as  shown  in  Figure  3,  in  the  spherical 
coordinate  system  commonly  used  for  antennas.  For  an  antenna  at  the  ■•igh  , 
0  nil  <1>  to  point  P(X,Y,Z)  are 


(-)  =  tan 


10  0  <  7T  radians) 


(3  4) 


7 r  tan 


'  X\ 

lY/ 


(0  ^  'I*  ^  2  7 r  radians). 


(.1  m 
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FIGURE  3  LOOK-ANGLE  ©  AND  4>  IN  SPHERICAL  COORDINATE  SYSTEM 
FROM  AN  ANTENNA  AT  ORIGIN  TO  P  (X,  Y,  Z) 


The  gain  functions  are  approximated  by  a  two-level  mouel.  The  3  dB 
limits  of  the  antenna  main  beam  are  defined  by  the  angles  0],  0p,  <t>;, 
and  <l>  p  .  The  gain  function  is  then 


Gl©,^) 


gmsl 


for  ©i^©^  ©2 
and  <  (|>  <J>2 

otherwise 


(3  6) 


In  Equation  (3.6)  Gmb  is  the  gain  of  the  main-bean,  and  the  Gmsl  is  the 
mean  value  of  the  side  lobe  gains.  When  G^b  and  GjqgL  ire  unknown  they  are 
computed  as  follows:  (Ref  ?) 


and  (Ref  3) 


GMB  "  10  lo9l0 


477 

(©2  -  e1i(,i>2  - 


(3.7) 


^MSL  '  +  0-5  Gmb 


(3  8) 
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special  situation  exists  for  downward -pointing  antennas  in  that  @i, 
0  ,  <t>  1  and  $  'i  cannot  be  defined  as  above.  In  this  case,  now  main-beam 
angles,  0 and  0  ■>  f  are  defined  in  the  Y-Z  (pitch)  plane  and  <l>  j. ,  and 
<t>  in  the  X-Y  (roll)  plane  as  shown  in  Figure  k.  New  look  angles  are 
then  computed  using 


0  =  tan  1  [Da  sinBsin  4>/Y] 


(3.9) 


’=  tan  1  (Da  sinQcos  4>/Y) 


(3.10) 


where  DA  -  >/x  2  +  Y2  +  Z2 

The  gains  are  determined  using  the  new  angles  in  Equations  (3*6),  (3*7) >  and 

(3.9). 


Y 


FIGURE  4  LOOK-ANGLE  0  AND  <t>  ’  FOR  THE  SPECIAL  CASES  OF  A  DOWNWARD- 
POINTING  ANTENNA.  IS  THE  ANGLE  IN  THE  X-Y  (ROLL)  PLANE, 

AND  0  IS  THE  Z-Y  (PITCH)  PLANE 
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3- 2. 4. 2  Distance  Evaluation 


The  antenna  separation,  D  in  Equation  (3-2),  is  determined  using  the 
aircraft  model  shown  in  Figure  5«  Antenna  locations  are  specified  in  the 
usual  fuselage  station  (FS),  waterline  (WL)  and  butt-line  (BL)  system  used 
for  aircraft.  For  analysis,  these  locations  are  converted  by  the  computer 
into  cylindrical  coordinates  (p,  8,  Z)  ( p  and  Z  are  in  inches,  6  is  in 
radians).  The  Z-axis  of  the  cylindrical  coordinate  system  is  the  centroid 
of  the  aircraft  and  is  displaced  vertically  from  the  WL  =  0  line,  a  distance 
denoted  as  WLC.  The  base  of  the  conical  nose  section  is  specified  by  FSN. 


Y-Axis 


D  is  determined  by  using  the  combination  of  straight  lines,  conical 
spirals,  and/or  cylindrical  spirals  that  gives  the  shortest  distance 
between  the  two  antennas  over  the  aircraft  surface.  The  combination  used 
depends  on  the  relationship  of  the  two  antennas  to  the  fuselage  and  wings. 
With  both  antennas  on  the  fuselage,  a  conical  spiral  is  used  when  one  or 
both  antennas  are  forward  of  FSN,  and  a  cylindrical  spiral  is  used  when 
both  are  aft  of  FSN. 

For  the  conical  spiral,  the  distance  between  two  points, 

P1  (  pl>  6l>  zl)  and  pp  (  p2>  e2>  'L2>  is 
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D 


con 


0,c-02  d 


a2  +  b2 
2a 


In 


afl1  +  c 
ad2  +  d 


P2  -P, 

where  a  =  - 

(3.11) 


c  =ya2  +  b2+a2fl2 
d  J  a2  +  b2  +  a2  0* 


For  the  cylindrical  spiral  (  P\,  =  Pp  =  Pf) >  the  distance  between  two  point 
is 


D 


cyl 


■/ 


P  <0, 


+  <Z2-Z,^ 


(3.12) 


Straight  line  distances  are  found  using 

Dsl  "  p2  2p,  p2  cos  (02  -  0,)  +  (Z2  -  Z,)2  (3.13) 

./hen  one  antenna  is  off  the  fuselage  but  the  other  is  on,  as  illustrate! 
in  Figure  6,  D  is  a  combination  of  straight  line  and  cylindrical  spiral 


segments. 

The  tangent  point  on  the  fuselage  for  which  D  is  minimum 
(  ■  turning  po,  Pp,  Zp  is  the  point  off  the  fuselage) 

[  v4>I-  P?  1 


Pj  02  t  tan  1 


t  coordinate 


(3.14) 


(The  appropriate  sign  is  used  that  produces  minimum  D.)  The  Z  coordinate 
of  the  tangent  point  is 

Pf  \0 f  P,  I  Z2  +  Zf  \/p2  —  Pf 

Zf  - -  ■  -  -  (3.15) 

Pi  10,-0,  I+VP2-P? 
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now  becomes  p_  in  Equation  3*13  and  (0^,  z^)  become  (Q Z^)  respectively 
when  calculating  the  straight  line  segment. 


Antenna  2 
(^2'  ^2'  ^ 


Z  Axis 
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FIGURE  6  DISTANCE  BETWEEN  TWO  ANTENNAS  WHEN  ONE  IS  LOCATED  OFF 
FUSELAGE  AND  A  PORTION  OF  THE  PATH  IS  AROUND  THE  FUSELAGE 


.■/hen  both  antennas  are  off  the  fuselage  and  the  shortest  path  between 
them  is  around  the  fuselage,  D  consists  of  two  straight  line  segments  and 
one  cylindrical  spiral  segment  as  shown  in  Figure  7.  In  this  case,  there 
are  two  tangent  points.  The  6  coordinate  for  the  tangent  points,  0.^ 
and  0^2)  are  computed  using  Equation  (3*1^)  above,  and  the  Z  coordinates 
are 

B  (Z,+  Z2>  +  Z2  A, 

z  =  -  (3.16) 

B  (B  +  A1  +  A2) 


where  A,  =y -  pf 

A2  "Vp2~  Pf  (3.17) 

B  =  Pf  l$ti  —  $,2  I 

p(  =  radius  of  the  fuselage. 


Z,2  = 


Z ,  B  +  Z2A, 
B  +  A, 
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FIGURE  7  DISTANCE  BETWEEN  TWO  ANTENNAS  WHEN  BOTH  ARE  OFF 
THE  FUSELAGE  AND  THE  PATH  IS  AROUND  THE  FUSELAGE 


la  the  event  that  a  wing  intersects  the  propagation  path,  the  wing  must 
be  included  in  the  distance  calculation  as  well  as  the  shading  calculation. 

The  propagation  path  around  the  wing  edge  must  be  determined  such  that  its 
total  length,  D,  is  a  minimum.  First,  the  forward  or  aft  wing  edge  is 
selected,  depending  on  which  path  will  have  the  shortest  total  distance. 

This  is  determined  by  calculating  the  intersection  of  the  wing  and  a 
straight  line  between  the  two  antennas.  Once  the  wing  edge  is  determined, 
an  iterative  technique  is  used  to  find  the  point  on  that  edge  such  that 
the  propagation  path  is  n  minimum.  A  typical  example,  with  both  antennas 
off  the  fuselage  is  shown  in  Figure  8.  For  this,  the  distances  are 
determined  from  the  wing  point  to  each  antenna  using  the  techniques  described 
above.  The  result  is  the  minimum  distance  D. 

3 . 2 . h . 3  Fuselage  Shading  and  Diffraction 

Where  a  portion  of  the  propagation  path  is  around  the  fuselage,  or  any 
curved  surface,  allowance  must  be  made  for  shading  effects.  An  equation  for 
the  fuselage  shading,  as  a  nortion  of  the  shading  factor  SF  in  Equation  (3-2), 
was  derived  from  Hasser.Jiam  and  Ishimai'u  (Ref  4).  In  their  work,  a  function 
was  given  relating  propagation  around  an  infinite  conducting  cylinder  to 
that  over  a  flat  plane.  The  function  can  be  reasonably  approximated  by 
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-A 

<T?A  +  f ) 


SF 


C 


where  A  = 


V  = 

t  = 


5.478  x  10~3,for  A  <  26 
3.340  x  10_3,for  A  >  26 

0.5083,  for  A  <  26 
0.5621,  for  A  >  26 


(3.181 


and  SFC  • 

S'  I 

Xs  = 


o 


c 


fuselage  (cylindrical)  shading  factor  (dB) 
radius  of  cylinder  (meters) 

angle  around  cylinder  of  propagation  path  (radians) 
wavelength  (meters) 

distance  of  cylindrical  segment  of  propagation  path  (meters) 


The  error  for  this  approximation  Is  less  than  2  percent  In  the  normal  range 
of  iso  (A  <  100). 
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FIGURE  8  TYPICAL  PROPAGATION  PATH  WITH  WING  SHADING 
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3 . 2  J< . U  Wing  Shading  and  Diffraction 


When  a  portion  of  the  propagation  path  la  around  the  wing,  or  any  sur¬ 
face  edge,  allowance  must  be  made  for  diffraction  effects  (Ref  5  and  6).  An 
incident  electric  field  can  be  related  to  a  diffracted  electric  field  using 
optic  theory  arid  lowering  the  frequency  range  to  radio  frequency  (Ref  5). 
From  Figure  9,  the  diffraction  attenuation  function  for  the  wing  can  be 
calculated  using 


/ 

SFW  =  20  lofl1Q 


where 


(ox-ar)i  esc  (ax+ar) 
2  sin 


and  SFW  =  Wing  shading  factor  (dB) 

£x  =  I  ncident  electric  field  from  transmitter  at  wing  point 

=  Received  electric  field  at  receiver  antenna 

ax  =  Angle  between  Y'  axis  and  incident  wave  propagation  vector 

aT  =  Angle  between  Y'  axis  and  receptor  wave  propagation  vector 

0  =  Azimuthal  angle  of  incident  wave  propagation  vector  (in  wing  plane) 

k  =  2irf/c  =  Propagation  constant 

c  =  Velocity  of  Light  (3x10®  m/sec) 

D  =  Distance  from  wing  point  to  receptor  point  (meters) 
wr 


(3.19) 


The  total  shading  factor  in  Equation  (3*2)  is  then  the  sura  of  the  edge 
and  cylindrical  surface  factors. 

SF  =  SFC  +  SFW  <3  20> 

3.3  analytic  assumptions  and  approximations 

The  mathematical  model  includes  assumptions  and  approximations.  This  is 
necessary  when  modeling  complex  or  unknown  functions,  such  as  antenna  pattern 
for  computer  analysis.  These  can  be  divided  into  receiver,  transmitter, 
antenna,  and  default  categories. 

3.3.1  Receivers 

The  following  assumptions  and  approximations  are  included  in  the 
receiver  sensitivity  model: 

o  The  receiver  sensitivity  is  constant  throughout  its  operating 
frequency  range  (minimum  discernible  signal). 
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o  The  selectivity  skirts  are  log-linear. 

o  There  are  no  receiver  responses  below  Far  and  above  F^r  (no  spurious 
responses;  significant  image  frequency  responses  are  inputted  as 
separate  receiver  responses). 

o  There  is  no  signal  discrimination;  that  is,  any  signal  exceeding 
the  threshold  level  will  appear  at  the  receiver  output. 


Y' 
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FIGURE  9  DIFFRACTION  COORDINATE  SYSTEM  WHEN  THE  PROPAGATION 
PATH  IS  AROUND  THE  WING  EDGE.  ATTENUATION  FUNCTION  IS 
CALCULATED  USING  ANGLES  IN  THIS  COORDINATE  SYSTEM 


3. 3- 2  Transmitters 

The  following  assumptions  and  approximations  are  included  in  the 
transmitter  spectrum  signature  model: 

o  The  power  output  is  constant  throughout  the  fundamental  operating 
frequency  range  (i.e.,  transmitter  rated  peak  power). 

o  The  skirts  are  log-linear. 

o  The  only  spurious  outputs  are  harmonics. 

o  The  spectrum  shape  of  each  harmonic  is  the  same  as  that  of  the 
fundamental. 
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3.3*3  Antennas 


To  approximate  antenna  patterns,  the  ECAC  two  level  model  was  used. 

This  model  assumes: 

o  The  gain  in  the  main  beam  (G^)  is  constant  for  all  angles  within 
the  3  dE  points. 

o  The  gain  outside  of  main  beam  is  constant  and  equal  to  Gj^l  (mean  of 
sidelobe  gain). 

o  The  relationship  of  the  approximation  model  to  a  typical  radar 
antenna  pattern  is  shown  in  Figure  10.  Three-dimensional 
approximation  patterns  for  typical  radar  and  communication  antennas 
are  illustrated  in  Figures  11  and  12,  respectively. 

o  The  antennas  are  100  percent  efficient  and  there  is  a  perfect 
impedance  match  to  the  load  or  source. 


0° 


FIGURE  10  TWO-LEVEL  APPROXIMATION  MODEL  FOR  DIRECTIONAL  ANTENNAS 
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FIGURE  11  ANTENNA  MODEL  APPLIED  TO  TYPICAL  RADAR  ANTENNA 


FIGURE  12  ANTENNA  MODEL  APPLIED  TO  TYPICAL 
COMMUNICATIONS  ANTENNA  WITH  DOUGHNUT  PATTERN 
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3.3*4  Default  Assumptions 


The  above  assumptions  and  approximations  apply  to  the  basic  mathematical 
model  and  are  used  in  all  calculations.  Also  included  in  the  procedure  are 
default  assumptions  for  use  only  when  certain  equipment  parameters  are  not 
known  or  not  specified.  These  are  as  follows: 

o  Transmitter  harmonics  are  60  dB  below  the  fundamental,  and  the 
highest  significant  harmonic  is  the  third.  (The  program  may 
accommodate  to  the  sixth.) 

o  Receiver  selectivity  skirts  conform  with  Figure  23,  page  26,  of 
MIL-I-6181D  (to  the  -80  dB  point). 

o  Transmitter  spectrum  skirts  are  as  defined  above. 

o  The  antenna  main  beam  gain  is  calculated  using  Equation  (3*7)  • 

o  The  antenna  mean  side  lobe  gain  is  calculated  using  Equation  (3*8)* 

3.4  COMPUTER  PROGRAM  DESCRIPTION 

3.4,1  Input  Data  Required 

This  is  a  brief  summary  of  the  input  parameters  required  by  ATACAP. 

The  information  is  categorized  into  aircraft  and  subsystem  data.  Subsystem 
data  is  further  divided  into  transmitter,  receiver,  and  antenna  data. 

Aircraft  Data.  This  data  describes  the  physical  dimensions  of  the 
aircraft  and  fits  the  cylindrical  model  to  the  actual  vehicle  shape.  There 
is  also  an  option  for  using  a  round  or  flat  bottomed  cylinder.  All  of  the 
data  is  in  inches  and  relative  to  butt-line  (BL),  water-line  (WL),  and 
fuselage  station  (FS)  as  commonly  used  on  aircraft.  The  aircraft  data 
required  are  summarized  as  follows  (see  Figure  5): 

o  Centroid  WL 

o  Conical  Nose  Limit  F5 

o  Bottom  WL  (for  flat-bottomed  cylinder) 

o  Location  of  Wing  Corners  (WL,  BL,  FS) 

Subsystem  Data.  For  the  purpose  of  this  analysis,  a  subsystem  is  defined 
as  a  group  of  single  band  transmitters  and  receivers  using  a  series  of  common 
antennas .  (Multiple  bands  and  spurious  responses  are  specified  as  separate 
receivers  and  transmitters.)  The  following  data  is  used: 

o  Transmitter  Data 

o  Frequency  range  covered 
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o  Skirt  frequencies 
o  Output  power 
o  Harmonics  (up  to  sixth) 
o  Antennas  used 
o  Emission  bandwidth 
o  Receiver  Data 


o  Frequency  range  covered 
o  Skirt  frequencies 
o  Sensitivity 
o  Antennas  used 
o  Bandwidth 
o  Antenna  Data 


o  Antenna  type 

o  Location  on  aircraft  (BL,  WL,  FS) 
o  Main  beam  gain 
o  Mean  side  lobe  gain 
o  Main  beam  limit  angles 
3.4.2  Analysis  Procedure  (Algorithm) 

The  procedure  used  in  the  program  to  implement  the  mathematical  model 
is  illustrated  in  the  flowchart  of  Figure  13*  The  flowchart  is  more  easily 
understood  if  the  two  analysis  steps  are  kept  in  mind: 

o  Frequency  coincidence 

o  EMI  margin 

The  rest  merely  involves  data  routing  and  manipulation. 

First,  all  the  data  for  one  aircraft  is  read  from  punched  cards.  The 
computer  examines  the  data,  decodes  the  alphanumeric  codes,  and  inserts  the 
default  values  where  needed.  Also,  the  hipest  security  classification  of 
the  subsystems  is  determined  so  that  it  can  be  printed  at  the  top  and  bottom 
of  each  page.  A  summary  of  the  input  and  default  data  is  then  printed. 
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FIGURE  13  (CONCLUDED)  ATACAP  ALGORITHM  FLOWCHART 
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The  analysis  involves  selecting  each  receptor  and  determining  what 
sources  will  interfere  with  it.  A  transmitter  and  receiver  are  selected 
from  the  data.  Their  frequency  ranges,  including  the  significant  harmonics 
of  the  transmitter,  are  compared  for  coincidence.  If  no  coincidence  is 

found,  compatibility  is  assumed,  and  they  are  not  analyzed  further.  If  there 

is  coincidence,  the  EMI  margin  is  determined  by  calculating  the  power  received 
at  the  receiver  front-end  and  comparing  it  to  the  receiver  threshold.  Worst  - 
case  EMI  margins  are  calculated  for  all  antenna  combinations  using  the 
shortest  path  along  the  aircraft  surface.  All  cases  with  margins 

exceeding  the  value  specified  in  the  input  data  are  listed  in  the  computer 

output.  This  procedure  is  repeated  for  all  subsystem  combinations  on  the 
vehicle.  The  EMI  margin  calculations  include  shading  due  to  the 
aircraft  surfaces,  diffraction  around  wings,  free  space  transmission  losses, 
and  coaxial  cable  losses. 

The  first  time  an  antenna  pair  is  encountered,  the  relevant  propagation 
data  such  as  distances,  angles,  antenna  gains,  etc.,  is  stored.  When  the 
same  antenna  pair  is  again  encountered  at  a  different  frequency  (harmonic) , 
only  the  frequency -dependent  portions  of  the  propagation  equations  are 
calculated.  This  saves  much  redundant  calculation  and  speeds  execution 
time  considerably. 

ATACAP  contains  a  "new  systems"  feature  which, for  reasons  of  clarity, 
was  omitted  from  the  flowchart.  This  feature  is  used  when  a  set  of  new 
subsystems  is  added  to  an  existing  vehicle.  Since  the  existing  avionics  are 
assumed  compatible  with  each  other,  they  are  not  compared  for  EMI,  but  the 
new  subsystems  are  compared  with  each  other  and  with  the  old  systems. 

3.5  PRELIMINARY  VERIFICATION 

The  ability  of  ATACAP  to  successfully  predict  EMI  was  checked  by 
comparing  early  predictions  with  later  tests  on  F-4  aircraft  containing 
complex  electronic  subsystems  and  numerous  antennas  including  ECM/RKAW 
equipment.  Examination  of  Air  Force,  Navy,  and  MCAIR  EMC  test  records 
showed  112  tests  were  conducted  against  which  the  ATACAP  predictions  could 
be  compared.  Of  the  interference  cases  found  in  the  tests,  only  2  percent 
of  these  cases  were  not  predicted,  and  only  18  percent  of  the  predicted 
interference  cases  did  not  appear. 

The  numerical  accuracy  of  the  model  has  also  been  checked  by  comparing 
the  ATACAP  predicted  antenna  coupling  with  132  measurements  on  a  KC-135 
aircraft  made  by  the  Air  Force  Avionics  Laboratory.  The  results  are  summarized 
in  a  histogram  in  Figure  l4. 

3.6  APPLICATION  EXAMPLE  AND  SAMPLE  OUTPUT 

A  typical  use  of  ATACAP  was  a  study  performed  for  the  U3AF  (FSCP  10R2 
M0D23b2) .  In  this  study,  a  number  of  transmitting  and  receiving  subsystems 
were  evaluated  for  compatibility  on  the  F-4c,  F-4d,  and  F-4e.  A  typical 
configuration  is  shown  for  the  F-4-D  in  Figure  15. 
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The  data  for  the  subsystems  of  Figure  15  was  fed  into  ATACAP,  and  a 
sample  of  the  computer  output  is  shown  in  Figure  16.  The  source  and^ 
receptor  subsystems  are  identified  in  the  output  for  each  predicted  hMI 
situation.  The  EMI  margin  is  shown,  and  other  data  is  presented  to  aid 
the  user  in  evaluating  the  interference  and  determining  how  to  prevent  Lt. 
The  maximum  EMI  frequency  is  shown  with  the  minimum  and  maximum  coincidence 
frequencies.  Also  given  are  the  free  space  propagation  1oss,  the  fuselage 
shading,  and  the  wing  shading,  all  in  dB.  Codes  arc  given  which  show  the 
type  of  frequency  coincidence  and  the  type  of  propagation  path  between  the 
antennas.  Finally,  an  output  is  provided  for  convenience  to  show  whether 
or  not  both  the  receiver  and  transmitter  have  inter-subsystem  blanking 
provisions.  The  results  of  this  example  analysis  are  summarized  in 
Table  I. 


TABLE!  ANALYSIS  RESULTS 
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1  Numbers  are  EMI  margin  in  dB  (amount  by  which  receiver  threshold  is  exceeded) 


2  <t>  shows  EMI  margin  less  than  —20  dB 
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F-4E 


ATACAP  COMPUTER  PROGRAM  TEST  RUN 


RECEPTOR  --  SUBSYSTEM  =IMAGPOO  ,  RCVR  =LO  BNO 


SOURCE  XMTR  HARM-  HA-X  -E-MI-  MIN  GOIMC - H AX-00 IHO - COINCIDENCE  8LKG 

SUBSYSTEM  IOENT  FREQ-  MHZ  FREQ-  MHZ  FREQ-  MHZ  XMTR  RCVR 


UHF  COMM  COMXTR  3  IjOQ.OOO  986.639  1362.271  MAIN  MAIN  NO 


EMI  RCVR  ANTENNA  XMTR  ANTENNA  *— PROPAGATION--*  PROP  CODE 
HARO  IOEHT  tOOE - IOEN-T - UOHE - FR  -SP - FOSE - HI  HO - R-H  -XM  HR 

-36.5  INPVOR  MAIN  NGSSOR  -MAIN  -42.  5 - 0^0 - 0.0 - 9-11— 0  -0- 

-34.3  INPYOR  MAIN  AOF  MAIN  -43.3  0.0  0.0  9000 

--3T.5 - Q  SPY  OR — MAIN - NOSEQR - MAIN - -4  6-.  6 - CU43 - (UO - 9-0 - 0 - 0— 

-37.1  OBPYOR  MAIN  AOF  MAIN  -46.1  0.0  0.0  9000 

-35.2  CENTOR  MAIN  AOF  MAIN  -44.2  0.0  0.0  900  0 


—  -SOURCE- - XHT-R - HARM - 4AX-SHI - MIN  GOIN-C - MAX-  CRtNC — COINCIDENCE  64 KG 

SUBSYSTEM  IOENT  F REQ-  MHZ  FREQ-  MHZ  FREQ-  MHZ  XMTR  RCVR 


TRANSPOR  TSPXTR  1  1.90.003  942.294  1237.706  MAIN  MAIN  YES 


EMI  RCVR  ANTENNA  XMTR  ANTENNA  •— PROPAGATION — *  PROP  CODE 
-  — ■  - MARC - 1  GENT - LORE - TRENT - CORE - FR  SP F-USfc HENC P  H  XH  -NR - 

57.9  INPYOR  MAIN  - TGPOGK  MAIN  -40.3 - 7.6 — 46.0 - 7  4  21  43 - 

61.3  OBPYOR  MAIN  TOPBCK  MAIN  -43.8  -3.4  -12.3  7  1  21  40 

52.1  CENTOR  MAIN-  TOPRCX - MAIN - -46.4  -3.3  -18.9 - 7  1  El  21 - 


SOURCE  XMTR  HARM  MAX  EMI  MIN  COINC  MAX  COINC  COINCIDENCE  BL KG 

SUBSYSTEM  IOENT  FREQ-  MHZ  FREQ-  MHZ  FREQ-  MHZ - XHTR - RCVR - 

TACAN  TACXTR  1  1,25.000  686.102  1263.126  ■  MAIN  MAIN  YEC— 
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FIGURE  16  SAMPLE  ATACAP  COMPUTER  OUTPUT  FORMAT 
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SFCTION  b 


WIRF-TO-WIRF  COMPATIBILITY  ANALYSIS  PROGRAM  (WTWCAP ) 


h . 1  BACKGROUND 


Achieving  wire-to-wire  compatibility  on  modern,  high-performance  air 
vehicles  is  a  formidable  task.  An  F-k  aircraft,  for  example,  contains  over 
16,000  wire  segments.  In  advanced  vehicles,  more  wires  are  needed,  and  higher 
density  bundles  with  smaller  wires  and  thinner  insulation  are  being  used.  An 
analysis  tool  is  needed  that  can  rapidly  and  reliably  indicate  the  FMI  in 
such  bundles. 

The  current  state-of-the-art  was  determined  by  reading  all  of  the  avail¬ 
able  reports  (Refe  7  -  14).  Investigation  was  conducted  into  applicable 
electromagnetic  field  theory,  circuit  theory,  transmission  line  theory,  and 
numerical  analysis  methods.  From  this,  mathematical  models  and  analysis 
techniques  were  tested  on  the  computer  and  compared  to  measured  data. 

h.2  PFTAILED  DESCRIPTION  OF  WTWCAP  PREDICTION  MODEL 

t.2.1  Basic  Model 


Coupling  between  wires  is  a  multi-wire  transmission  line  problem  because 
the  coupling  mechanisms  are  distributed  along  the  wires.  In  general,  RF 
magnetic  and  electric  field  effects  are  interrelated  and  cannot  be  separated. 
However,  at  frequencies  where  the  wavelength  Is  large  compared  to  the  bundle 
lengths,  a  lumped-parameter  circuit  quite  accurately  approximates  the  dis¬ 
tributed  parameter  circuit.  Such  a  lumped  parameter  pi  circuit  is  illustrated 
in  Figure  17  for  two  wires  above  a  ground  plane. 


FIGURE  17  TWO  WIRE  COUPLING  WITH  SINGLE  PI  LUMPED  PARAMETER  CIRCUIT  MODEL 
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As  the  frequency  is  increased,  accuracy  is  lost  because  the  lumped 
parameters  no  longer  represent  the  distributed  parameters.  However,  if  the 
number  of  circuit  pi  sections  is  increased,  as  shown  in  Figure  18,  accuracy 
can  be  maintained  to  as  high  a  frequency  as  desired.  This  can  be  seen  in 
Figure  19  which  shows  the  transfer  function  for  one,  two,  and  three  pi 
sections.  A  mathematical  Justification  for  this  is  presented  in  the  Appendix. 


Rw  L  Rw  L 


FIGURE  18  TWO  PI  CIRCUIT  MODEL 


100  K  1  M  10  M  100  M  1  G 


Frequency  •  Hertz  GP71  0386  22 

FIGURE  19  TRANSFER  FUNCTION  FOR  1,2,  and  3  PI  LUMPED  PARAMETER 
CIRCUIT  MODELS  (10  FT  WIRES) 


30 


Frequency  domain  is  used  for  the  analysis  because  it  allows  the  use  of 
a  varying  number  of  pi  sections.  A  single  pi  circuit  is  used  at  low  fre¬ 
quencies,  a  double  pi  circuit  at  higher  frequencies,  and  so  on.  Thus, 
maximum  computer  efficiency  is  achieved  since  the  least  complicated  circuit 
is  used  consistent  with  accuracy.  Also,  frequency-dependent  wire  parameters, 
such  as  skin  effect  on  wire  resistance,  can  be  included.  Furthermore,  the 
complex  matrix  equations  resulting  from  frequency  domain  analysis  can  be 
more  readily  implemented  on  the  computer  than  the  numerical  methods  required 
to  solve  simultaneous  differential  equations  in  the  time  domain.  Finally, 
frequency  domain  is  directly  relatable  to  EMC  tests  and  specifications  so 
that  they  can  be  applied  directly  to  the  analysis.  If  time  domain  parameters 
of  the  interference  are  required,  they  can  be  obtained  from  the  complex 
numbers  by  the  inverse  Fourier  transform. 

For  the  wire-to-wire  problem,  simultaneous  analysis  is  used  since  it  is 
more  efficient  on  the  computer  and  can  account  for  second  and  higher  order 
effects.  (Higher  order  effects  result  from  signals  coupling  onto  one  wire, 
these  signals  then  coupling  onto  other  wires,  these  signals  coupling  onto  still 
more  wires,  and  so  on.  In  other  words,  every  wire  effects  every  other  wire.) 
By  solving  the  entire  circuit  simultaneously,  these  effects  are  automatically 
included.  Our  investigation  showed  it  is  practical  to  analyze  a  two  pi 
circuit  of  60  conductors  on  the  CDC  6600  computer  (250Kg  memory).  Larger 
bundles  can  be  divided  and  run  in  segments,  or  only  the  worst-case  sources 
and  receptors  can  be  analyzed. 

k.2.2  Circuit  Parameters 


The  parameters  involved  (Figure  17)  are  the  inter-wire  capacitances  and 
mutual  inductances,  wire  self-inductance,  wire  resistance,  and  wire  capaci¬ 
tance  to  ground.  These  are  divided  into  frequency-independent  and  frequency- 
dependent  categories.  In  the  following,  all  units  are  in  the  rationalized 
MKS  system. 

U.2.2.1  Constant  Parameters 

The  frequency- independent  parameters  include  the  external  self-induc¬ 
tances,  mutual  inductances,  capacitance  to  ground,  and  inter-wire  capaci¬ 
tances  . 

Self- inductance  of  wires  and  shields  are  best  determined  by  considering 
the  effects  of  internal  and  external  magnetic  fields.  Skin  effect  involves 
only  the  internal  fields,  so  that  the  external  fields  are  frequency  inde¬ 
pendent.  The  external  self-inductance  of  a  wire  of  radius  r.  ,  length  1, 
height  h,  carrying  a  current  L^,  as  shown  in  Figure  20,  is  defined  by 
(Ref.  12) 
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is  the  magnetic  flux  density  resulting  from  1^,  and  A^  is  the  surface 
between  the  wire  and  its  image  in  the  ground  plane.  Evaluating  Equation  (4.1) 
for  Figure  20  leads  to 
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2  IT 
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(4.2) 


henries/meter  where  /x  is  the  permeability  of  the  medium.  Since  Equation  (4.2) 
only  involves  the  external  fields,  it  also  applies  if  the  wire  is  hollow. 

Thus,  this  applies  both  to  solid  wires  and  to  shields. 


The  situation  for  the  mutual  inductance  between  two  wires  is  shown  in 
Figure  21.  Mutual  inductance  is  defined  by  (Ref.  12) 
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1 2  " 


(4.3) 


where  is  the  magnetic  flux  density  from  wire  1  and  A2  is  formed  by  wire  2 
and  its  ground  plane  image.  Note  that  reciprocity  applies  so  that  M-^  = 

For  the  mutual  inductance  to  a  wire  or  shield  to  another  wire  or  shield. 
Equation  (4.3)  results  in 


m,2 


4  7T 


(4.4) 


FIGURE  20  EXTERNAL  SELF  INDUCTANCE  OF  WIRE  OR  SHIELD 
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FIGURE  21  MUTUAL  INDUCTANCE  BETWEEN  WIRES  OR  SHIELDS 


hen::ies/meter  where  s  is  the  separation  between  centers.  For  the  case  of  a 
wire  to  its  own  shield,  as  shown  in  Figure  22,  Equation  (4.3)  becomes 

H  2h 

Mws  =  ~  lo9e  — 

4  7r  r0 

where  r0  is  the  outer  radius  of  the  shield. 

Capacitance  parameters  are  determined  from  the  basic  electrostatic 
charge  and  potential  relations.  For  the  case  of  two  wires  above  ground, 
as  shown  in  Figure  23,  the  relations  between  charge  and  potential  are 
(Kef.  12) 


01  ~  Pn  Q1  +  P12  ^2 
0 2=  P21  ^1  +  P22  ^2’ 


(4.6) 
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FIGURE  22  MUTUAL  INDUCTANCE  FROM  INTERNAL  WIRE  TO  SHIELD 


FIGURf  23  SITUATION  FOR  OPEN  WIRE  CAPACITANCE  CALCULATION 
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Here,  and  </>2  are  the  potentials,  and  Q,  and  Q2  are  the  charges  on  wires  1 
and  2,  respectively.  For  given  charges,  the  potentials  can  be  determined 
using  the  method  of  equipotential  surfaces  (Refs.  12,  13,  lM*  Referring 
to  Figure  23*  the  coefficients  of  Equation  (U»6)  are  as  follows: 

pn  =  cosh-1  <h  1  /r ■, >  (4  7) 

p22  =  cosh-1  (h2/r2)  (4-8) 


Pi2  -  P21 


s2  +  (rf  -  ri) 


2  sr. 


+  cosh  1 


(f?-r2) 


2  sr2 


cosh-1 


+  (rf-r2) 


2  s'r. 


-  cosh  - 1 


<'?■ 


rj) 


2  s'r^ 


(4.9) 


where  s'  is  the  diagonal  separation  to  the  image  and  is  given  by 

s'  =>/s2+(2h1)2-4sh  cos  6. 


(4.10) 


Then,  the  capacitance  relations  are 


C 


1  G 


27T6  pn 


det 


C 


2G 


lire  p22 


det 


p 

where  det  =  p-j^  p22  -  P12 
The  effective  permittivity. 


e 


eff 


'12 


-2TCeff  P12 
det 


is  the  determinate  of  the  p  matrix. 
«eff,  is  given  by  (Ref.  7) 


r  r,  +  r  " 

2 

1  2 

-1 

r.  +  r_ 

1  2  j 

1 

'2s  +  R,  +  R2‘ 

2 

-1 

2 

r,  +  r 

1 

2 

(4.11) 


(4.12) 


(4.13) 


(4.14) 


where  h^and  R2  are  the  wire  radii  including  the  insulation. 

In  the  above,  and  C2q  are  the  capacitances  to  ground,  and  CN2  is  the 
inter-wire  capacitance  in  farads/meter.  These  equations  apply  to  hollow  as 
well  as  solid  wires  and  therefore  can  be  used  for  shield-to-shield,  shield- 
to-wire,  and  shield-to-ground  capacitances. 

A  similar  approach  can  be  used  for  wires  within  a  shield.  The  shield 
is  used  for  reference  in  place  of  the  ground  plane,  and  capacitance  from 
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wires  within  the  shield  to  external  wires  and  to  the  ground  plane  are  assumed 
zero.  The  model  can  handle  up  to  4  shielded  wires.  For  a  single  shielded 
wire,  the  wire  to  shield  capacitance  is  (Ref.  12) 


2me 


farads/meter  where  rgi  is  the  internal  shield  radius  and  rw  is  the  inner 
conductor  radius. 

For  shielded  groups  of  from  two  to  four  wires,  the  following  procedure 
is  used:  The  potential  coefficient  matrix  terms  can  be  calculated  using 

I  cosh-1  (h2/rw)  for  i  -  j  =  0 


P.1  =  1 

loge  (D,  /D2 )  for  |  i  -  j 

odd 

(4  16) 

2a  +  2b  +  c 

i  log.  -  for  |  i  -  | 

2b  +  c 

I  even 

where 

b  =  x  <rw  + 

(4.17) 

In  Equation  (4.17),  x  is  given  in  Table  F"  -  ia 
the  wire  insulation  thickness.  Next, 

the  wire  radius,  tw  is 

hi  =  +  r  *  -  b2)/(2b) 

(4.18) 

where  H  .  is  the  shield  internal  radius.  Then, 
si 

a=\/h? 

(4.19) 

c  =  h,  -  b  -  a 

(4.20) 

h2  =  a  +  c 

(4.21) 

The  expressions  for  and  are  given  in  Table  II, 


TABLEH  SHIELDED  WIRE  CAPACITANCE  PARAMETERS 


No.  of 

Wires 

in  Shield 

X 

(Eq  4.17) 

D1 

d2 

2 

1.0000 

b  +  h1  +  a 

2b  +  c 

3 

1  1547 

v/(h1  +  a)2  +  b2  +  (h1  +  a)  b 

\J b2  +  (b  +  c)2  +  b  (b  4'  c) 

4 

1  4142 

/(h,  +  a)2  +  b2 

>/b2+  (b  +  c)2 

GP71  0386  19 
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(4.22) 


The  inverse  p  matrix  is  calculaced  such  that 

[cj  =  [pp1 

Then,  the  capacitances  are  found  using 

n 

C,s  =  2^S  (4  23) 

K=  1  ' 

C  =  -2mec  for  i  t  j  (4.24) 

ij  ij  1 

and  C.jj  is  the 


wires  and  shields, 
to  obtain  complex 


(4.25) 


where  the  primes  denote  differentiation  cf  the  Bessel  functions  with  respect 
to  their  arguments.  Breaking  this  into  real  and  imaginary  parts  leads  to 
expressions  for  the  resistance  and  internal  self-inductance  per  meter  of  wire. 
In  practice,  computations  of  these  involved  expressions  can  be  simplified  by 
fitting  them  with  approximation  curves.  In  the  computer  program,  these  fitted 
curves  are  much  faster  and  give  results  within  2  percent  of  the  above  expres¬ 
sion. 


CjS  is  the  capacitance  between  the  i™1  wire  and  the  shield, 
capacitance  between  the  i^*1  and  wires  in  farads/meter. 

b.2.2.2  Frequency-Dependent  Parameters 

These  parameters  involve  the  internal  impedance  of  the 

The  Helmholtz  equation  can  be  applied  to  this  problem, 
internal  impedance  of  a  solid  wire  (Ref.  12) 


Frequency-dependent  shield  parameters  can  be  analyzed  in  the  same 
fashion.  For  large  radius  of  curvature  and  large  shield  thickness  compared 
to  skin  depth,  the  following  simplified  results  can  be  applied  for  the  resis¬ 
tance  and  inductance  of  a  shield  (Ref.  12) 


where 

5  Uy'nifia 
1 


ts  shield  thickness 


and 


Ds  shield  inside  diameter 


(4.26) 


(4.27) 
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*4.2.3  Application  of  Model  to  Twisting  and  Shielding 


Twisting  and  shielding  are  important  methods  for  controlling  wire-to- 
wire  FMI.  Consequently,  the  program  must  have  an  accurate  method  for 
accounting  for  them. 

Since  most  shields  used  in  aerospace  vehicles  are  low  permeability, 
their  primary  effect  is  to  reduce  capacitive  coupling.  In  the  model,  there¬ 
fore,  capacitance  from  shielded  wires  to  external  wires  is  set  to  zero.  The 
low  permeability,  however,  does  not  prevent  magnetic  coupling,  and  the  mutual 
inductance  is  retained.  Shield  pigtails  are  modeled  as  inductances  of  0.0? 
/xhy  if  the  exact  impedance  is  not  known.  An  application  of  this  to  the  WTWCAP 
model  is  illustrated  in  Figure  2k. 


To  reduce  magnetically  coupled  EMI,  twisted  wires  are  used.  This  situ¬ 
ation  is  illustrated  in  Figure  25  in  which  currents  in  the  source  wire  I5j 
induce  equal  noise  currents  ( In )  in  the  two  receptor  wires.  In  practice, 
however,  the  receptor  wires  are  not  perfectly  balanced  with  respect  to  the 
source  wire.  For  this  reason,  the  mutual  inductance  in  the  model  from  t.hr 
source  wire  to  the  return  wire  is  reduced  slightly. 
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4.2.4  Wire  Separation  Considerations 

An  extensive  investigation  of  wire  positioning  within  a  bundle  showed 
that  the  exact  location  of  a  given  wire  cannot  be  determined  exactly  but  can 
only  be  estimated  on  a  probability  basis.  For  this  reason  and  since  worst- 
case  EMI  is  desired,  the  wire  parameters  are  calculated  assuming  minimum 
separation.  The  only  exception  is  for  twisted  groups,  as  previously  dis¬ 
cussed. 

4.3  COMPUTER  PROGRAM  DESCRIPTION 


The  Input  data  is  easily  obtained  from  information  usually  available 
to  the  user.  In  some  cases,  options  are  available  to  allow  the  user  to 
specify  data  in  the  form  in  which  it  is  available. 

Wire  Characteristics  Table  -  To  free  the  user  from  having  to  specify 
and  remember  nominal  wire  characteristics  for  each  type  of  wire,  a  table  of 
these  parameters  is  first  specified.  Each  type  of  wire  is  identified  by  a 
user  supplied,  alphanumeric  name  (e.g.,  22TFE  to  identify  22  gage,  teflon 


insulated  wires).  Wire  diameter,  conductivity,  insulation  thickness  and 
dielectric  constant,  and  shield  data,  where  applicable,  are  specified.  In 
specifying  the  bundle  data,  only  the  wire  name  is  used  to  identify  the  type, 
and  the  computer  obtains  the  nominal  characteristics  from  the  table. 

General  Bundle  Data  -  This  data  aids  in  identifying  the  bundle  in  a 
run  and  gives  general  parameters.  The  following  are  included: 

0  Bundle  designation  or  nomenclature 

°  Subsystem  and  unit  nomenclature 

0  Connector  nomenclature 

°  Bundle  length 

0  Heigi  t  above  ground 

Individual  Wire  Data  -  This  data  specifies  parameters  associated  with 
each  wire,  twisted  wire  groupings,  and  grounding.  The  following  are  included: 

o  Wire  type  designation 

0  Twisted  group  information 

°  Pin  designations 

0  Signal  reference  or  return 

°  Function  of  loads  (voltage  source,  voltage  sensitive  load,  current 
source,  current  sensitive  load,  grounded  wire,  etc.) 

0  Termination  impedances 

0  Shield  termination  data 

0  Source  and/or  susceptibility  spectra 

Spectra  data  -  Source  and/or  susceptibility  spectra  data  can  be  speci¬ 
fied  three  ways: 

0  Read-in  from  measured  data 

0  Built-in  EMC  specifications 

0  Built-in  spectrum  calculation  routines. 

If  the  spectrum  is  known,  such  as  from  measured  data,  it  can  be  specified  by 
giving  points  on  the  curve.  The  program  interpolates  (log-linearly )  between 
these  points.  If  EMC  military  specification  spectra  for  power  lines  is 
desired,  the  program  has  built-in  MIL-STD-704A  ripple  data  and  MIL-STD-W>1 
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data  for  conducted  susceptibility.  Finally,  if  commonly  encountered  time- 
domain  wave  forms  are  known,  the  program  will  calculate  the  spectrum  from 
the  wave  form  parameters.  These  waveforms  include: 

°  Step  function  (switch) 

°  Rectangular  pulse 

°  Trapezoidal  pulse 

°  Spike 

Additional  Inputs  -  Additional  inputs  provide  control  of  the  analysis 
and  output.  A  title  is  read  and  printed  in  the  output  for  run  identifica¬ 
tion.  A  printout  limit  is  specified,  and  cases  with  EMI  margins  below  that 
level  are  not  printed.  Details  of  a  given  frequency  range  and  a  plot  of  the 
calculated  EMI  can  be  ordered. 

4.3.2  Default  Assumptions  and  Options 


As  with  ATACAP,  many  parameters  that  are  required  are  unknown  or 
unavailable.  The  default  assumptions  and  options  contained  in  WTWCAP  are 
summarized  below. 

Unless  otherwise  specified,  the  following  bundle  parameters  are 
assumed: 

0  4  inch  height  above  ground 

0  10K  load  impedance 

0  0.02  fihy  impedance  from  the  shield  termination  to  ground 

If  the  bundle  length  is  unknown  and  the  box  locations  in  the  vehicle 
are  known,  the  length  is  approximated  by  calculating  twice  the  distance 
between  the  boxes. 

4.3.3  Analysis  Procedure  (Algorithm) 

An  algorithm  flowchart  is  presented  in  Figure  26.  The  input  data  is 
first  read  and  processed.  The  alphanumeric  codes  are  then  decoded  and  tested 
for  errors,  the  default  values  are  inserted,  and  the  appropriate  spectra 
data  are  calculated,  where  required.  The  input  data  summary  is  printed. 

The  constant  circuit  model  wire  parameters  are  then  calculated.  Where 
shields  are  encountered,  they  are  set  up  as  additional  wires  with  the 
appropriate  parameters.  For  shielded  wire  groups,  the  capacitance  to  ground 
and  capacitance  to  external  wires  are  set  to  zero.  Mutual  and  external 
self-inductances  are  then  computed  for  all  wires. 
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FIGURE  26  WTWCAP  ALGORITHM  FLOWCHART 
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FIGURE  26  (CONCLUDED)  WTWCAP  ALGORITHM  FLOWCHART 


The  EMI  analysis  is  initiated  by  setting  the  frequency  to  the  initial 
value  (usually  100  Hz).  Based  on  the  ratio  of  wavelength  to  bundle  length, 
the  number  of  circuit  pi  sections  is  selected.  The  circuit  model  package 
is  then  called  into  use. 

The  heart  of  the  program  is  the  circuit  model  package.  It  contains 
a  completely  general  scheme  for  setting  up  and  solving  by  nodel  analysis 
the  circuit  model  with  a  specified  number  of  pi  sections.  (The  number  of 
Pi  sections  is  limited  to  two  only  because  of  core  memory.  If  desired, 
more  pi  sections  could  be  run  with  fewer  wires  using  this  routine. ) 

First,  the  frequency-dependent  parameters,  internal  self-inductance 
and  wire  resistance,  are  calculated.  The  inductance  matrix  is  calculated 
and  inverted  for  use  in  the  nodal  admittance  matrix.  Next,  the  admittance 
matrix  is  calculated.  Included  are  the  wire  parameters  lumped  to  the 
specified  number  of  pi  sections.  The  termination  admittances  are  included, 
connected  to  ground  or  between  nodes  as  required. 

The  current  source  matrix  is  then  computed.  The  output  of  each  source 
at  the  present  frequency  is  determined  from  its  spectrum,  and  voltage  sources 
are  converted  to  current  sources.  The  appropriate  inter-node  and  ground 
connections  are  included. 

The  circuit  equations  are  solved  by  using  a  Gaussian  elimination 
routine  for  symmetrical  complex  matrices  to  obtain  the  complex  node  voltages. 
Then,  where  necessary,  inter-node  voltages  and  load  currents  for  current- 
sensitive  loads  are  calculated.  The  susceptibility  at  the  present  frequency 
is  determined  for  each  load,  and  the  EMI  margins  are  calculated  and  printed. 

The  circuit  analysis  is  repeated  for  all  desired  frequencies.  Plots, 
if  desired,  are  then  generated.  This  process  is  repeated  for  all  bundles. 

U .  3 .  U  Program  Outputs 

As  with  the  other  programs,  a  complete  summary  of  the  input  data 
including  any  default  values  is  first  printed.  The  constant  wire  parameter 
matrix  is  also  printed,  if  desired.  This  shows  how  the  shields,  shielded 
wire  groups,  etc.,  were  handled  by  the  program. 

Next,  the  frequency  by  frequency  analysis  is  printed.  For  each 
susceptible  receptor,  the  EMI  margin,  received  signal  level  (both  magnitude 
and  phase),  and  susceptibility  levels  are  printed  for  each  interference 
source.  The  complex  sum  of  EMI  from  all  sources  is  also  printed  in  each 
case . 


If  Plots  of  the  interference  are  desired,  the  computer  sorts  the 
calculated  data,  figures  the  scale  factors,  and  generates  semi-log  plots 
on  the  printer. 

If  they  were  requested,  detailed  frequency  range  analysis  outputs  are 
then  printed,  using  the  same  format  as  the  general  analysis.  As  of  this 
report,  no  plots  are  available  for  the  detailed  analyses. 
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U.4  PRELIMINARY  VER IFICATION 


A  few,  limited  tests  were  performed  to  verify  the  basic  assumptions 
used  In  the  analysis.  The  tested  bundle  is  shown  in  Figure  2J}  and  the 
measured  and  calculated  results  are  shown  in  Figure  28  for  the  twisted 
pair  from  the  28  volt  switch. 


FIGURE  27  PRELIMINARY  VERIFICATION  TEST  BUNDLE 


10K  100K  1 M  10M  100M 

Frequency  -  Hertz  GP71  0386  34 

FIGURE  28  MEASURED  vs  CALCULATED  EMI  INTO  TWISTED  PAIR 

(LOAD  END  OF  BUNDLE) 
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U.5  APPLICATION  EXAMPLE  AND  SAMPLE  OUTPUT 


A  typical  use  of  WTWCAP  can  be  demonstrated  by  the  case  shown  in 
Figure  29.  In  this  example,  a  digital  receiver  output  is  connected  to 
anon-board  computer  input.  The  minimum  spike  rejected  by  the  computer 
input  logic  circuitry  is  also  shown.  To  achieve  EMC,  the  following 
information  is  requii'ed; 

0  Is  it  feasible  to  use  this  bundle  in  the  aircraft  EM  environment? 
0  Is  the  given  spike  rejection  level  adequate? 

0  What  type  of  wiring  will  give  adequate  FMI  protection  with  mini¬ 
mum  weight? 

To  answer  these  questions,  a  bundle  was  devised  containing  the  worst-case 
sources  and  three  receptor  wire  configurations.  This  bundle  is  shown  in 
Figure  30.  In  this  way,  the  various  receptor  wires  could  be  tested. 


Digital  Circuit 


Minimum  Spike  Rejected  by  Logic 


FIGURE  29  WIRE-TO-WI RE  APPLICATION  EXAMPLE 
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Receiver  Aircraft  Computer 


The  results  of  the  WTWCAP  analysis  are  summarized  in  Figures  31  through 
35.  Figure  31  shows  the  induced  EMI  from  the  MIL-STD-TO^A,  115  volt,  hOO  Hz 
power  line.  As  shown,  the  shielded  pair  had  no  EMI,  the  twisted  pair  was 
marginal,  and  the  open  wire  has  definite  EMI.  Similarly,  the  SMI  from  the 
28  VDC  line  and  the  600  volt  relay  spike  are  shown  in  Figures  32  and  33, 
respectively.  A  portion  of  the  printed  computer  output  1b  shown  in  Figure 
3U,  and  a  sample  computer  plot  is  shown  in  Figure  35. 
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FIGURE  31  CALCULATED  INTERFERENCE  LEVELS  FROM 
MIL-STD-704A  115V AC  400  Hz  POWER  LINE 


FIGURE  32  CALCULATED  INTERFERENCE  LEVELS  FROM 
MIL-STD-704A  RIPPLE  ON  28VDC  POWER  LINE 
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FIGURE  34  WTWCAP  SAMPLE  NUMERICAL  OUTPUT  FOR  BUNDLE  OF  FIGURE  30 
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FIGURE  35  WTWCAP  SAMPLE  PLOT  OUTPUT  FOR  BUNDLE  OF  FIGURE  30 


SECTION  5 


ELECTROMAGNETIC  FIELD-TO-NIRE  COMPATIBILITY  ANALYSIS  PROGRAM  (FTWCAP) 
5.1  BACKGROUND 


The  coupling  from  external  electromagnetic  fields  onto  internal  wiring 
is  important  in  the  design  of  aerospace  vehicles.  Usually,  the  fields 
enter  the  vehicle  through  dielectric  apertures  in  the  vehicle  skin  and 
couple  onto  wires  immediately  adjacent.  These  apertures  include  radomes, 
canopies,  landing  gear  doors,  camera  windows  and  air  intakes.  FTWCAP 
presently  analyzes  only  coupling  onto  wires  from  on-board  transmitting 
antennas,  but  the  computation  method  is  completely  general.  It  can  be 
readily  adapted  to  analyze  EMI  from  other  RF  energy  sources,  such  as 
ground-based  transmitters,  lightning,  and  EMP. 

5.2  DETAILED  DESCRIPTION  OF  MODEL 

To  obtain  the  equivalent  plane -wave  RF  power  density  at  each  aperture 
from  the  on-board  radiating  antennas,  the  propagation  model  from  ATACAP  is 
used.  The  exposed  wires  are  assumed  to  be  adjacent  to  the  aperture,  and 
the  amount  of  RF  energy  coupled  depends  on  the  aperture  size  and  location. 

A  transmission  line  model  is  then  used  to  compute  the  currents  induced  in 
the  wire  loads  (Refs  15,  l6) .  Worst-case  electromagnetic  field  vector 
orientation  is  determined  and  used  for  the  calculation. 

5.2.1  Field -to-Transmission  Line  Equations 

Figure  36  shows  a  transmission  line  exposed  to  an  electromagnetic 
field. 


FIGURE  36  TRANSMISSION  LINE  ILLUMINATED 
BY  ELECTROMAGNETIC  FIELD 
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Integration  of  Faraday's  law  over  the  surface  of  the  transmission 
line  results  in  (Ref  13) 


f  ( VXE)  •  n  dS  =  -  —  fe  *  n 

wA,  at  J 


Defining  the  voltage  across  the  wires  at  a  given  z  as 


V  (z)  = 


leads  to  the  following  transmission  line  partial  differential  equations  with 
voltage  and  current  sources  included  (Ref  16) : 


3V(z)  f  b 

-  +ZI(z)  =  ja>  /  B'(x,z)dx 

9z  Jo 


3l(z)  f  b 

-  +YV(z)  =  jcjY|  Ex'(x,z)dx 

3z  Jo 


where  a  sinusoidal  time  dependence  is  assumed,  E*  and  are  the  incident 
fields,  and  Y  and  Z  are  the  transmission  line  admittance  and  impedance, 
respectively  (Ref  16) . 

Equations  (5*3)  and  (5.4)  can  then  be  used  to  derive  two  inhomogeneous 
second  order  differential  equations  for  the  transmission  line  currents 
and  voltages  (Ref  15) .  These  equations  are  solved  using  the  method  of 
variation  of  parameters  (Ref  17) . 

When  the  entire  transmission  line  is  not  exposed  to  the  incident 
field,  expressions  for  the  equivalent  impedance  of  the  unexposed 
transmission  line  and  the  physical  loads  replace  Zo  and  Zs  in  Figure  36 
and  the  equations.  The  resulting  current  in  Z0  is  then 


0  D  [Zc  cosh  (7^)  -  Zg  sinh  (76,] 


[J «,  ,E, 


{ b, z)  — E  (o,z)]  [Z  cosh  7  (C  -  z)  +  Zn  sinh  7  (  K  -  z ) ]  dz 

Z  C  Z  V.2  z 


+  (Zc  cosh  7(C2  -  C, )  +  Zg^  sinh  7  (C2  -  C,  )]J  E^  (x,  ^  )  dx 


"  Zc  I  Ex  ,x'  V  dx’ 
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Likewise,  the  current  in  Z$  is  given  by 


D  [Zr  cosh  7  (s-t2)  +  Z^  sinh  7  (s-  P,)] 


(/ 


(b,z)  -  (0,z)!  (Zc  cosh  7  (z  -  Cj)  +  Zg,  sinh  7  (z  -  fij)]  dz 


(5.6) 


IZc  cosh  (7 C2  )+  Zg  sinh  (7 fi2> )  /  E'Jx,  £2)  dx  +  Zc  /  E'x(x,  ^  )  dx 
’  ./  0  Jo 


where 

Zc  =  transmission  line  characteristic  impedance 

7  =  complex  propagation  constant 

b  =  separation  between  wires 

s  =  total  length  of  transmission  line 

^andfg  =  define  exposed  portion  of  line 

and 


D  =  {ZcZg,1  +  ZC2Zc)  cosh  7  <S2  -  ei )  +  (Zc  +ZC2  ZCi> sinh  7  -81 ) 


Zil  and 
edges.  They 


Zj£2  are  ‘the  termination  impedances  as  seen  at  the  aperture 
are  given  by 


(5.8) 


(5.9) 


These  equations  are  dependent  only  on  the  incident  electric  fields,  since 
the  scattered  fields  have  been  accounted  for  in  the  induced  transmission 
line  currents. 

5 .2.2  Multiple  Exposures 

The  above  analysis  can  be  easily  extended  to  wires  routed  close  to 
more  than  one  aperture.  Using  superposition,  the  analysis  is  used  to 
calculate  the  currents  from  each  exposure,  and  the  resulting  complex 
currents  are  summed. 


Zo  -  Zc  tanhTi, 

Zo.^Z,  - — 

Zc-Zotanh7«1 

Zs  +  Zc  tanh  7  (  S  -  C2) 

Zg2  =  Zc - 

Zc  +  z%  tanh7(s-e2) 
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5.2.3  Ground  Plane  Effects 


This  analysis  can  also  be  applied  to  circuits  with  a  ground  plane 
return.  In  this  case,  the  transmission  line  is  formed  by  using  the  ground 
plane  image  of  the  primary  wire.  The  above  equations  are  applied  using 
twice  the  height  above  ground  for  the  wire  separation. 

5.2.4  Shielding  Factors 

Allowance  must  be  made  in  the  model  for  the  overall  bundle  and 
individual  wire  shields.  Since  solid  shields  are  seldom  used  in  aerospace 
vehicles  because  of  weight  and  flexibility  considerations,  determining  the 
exact  degree  of  electromagnetic  field  attenuation  is  an  involved  task.  This 
determination  involves  the  size  and  shape  of  the  shield  weave,  the  number  of 
picks  or  openings  per  unit  length,  and  many  unknown  factors.  Therefore,  the 
shielding  factor  for  the  FTWCAP  model  was  based  on  empirical  data.  From 
many  measurements  on  all  types  of  braided  shields,  the  curve  shown  in 
Figure  37  was  selected  as  representative.  Two  typical  measured  curves 
for  a  single  overall  shield  are  also  shown.  The  shielding  effectiveness 
package  in  FTWCAP  is  separate  from  the  rest  of  the  program  and  can  easily 
be  updated  or  changed  if  other  shielding  effectiveness  is  desired. 
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FIGURE  37  MODEL  AND  TYPICAL  MEASURED  SHIELDING  EFFECTIVENESS 
CURVES  FOR  SINGLE  OVERALL  SHIELD 


5.2.5  Bandwidth  Considerations 


.  --  jrrTTrrt^f* 


Since  most  wire  load  susceptibilities  are  broad  oand,  the  received 
signal  must  be  in  terms  of  spectral  density  ( /u,V/ MHz)  before  they  can  be 
compared  and  an  EMI  margin  calculated.  Thus,  a  factor  is  used  to  convert 
the  relatively  narrow  band  transmitter  emission  to  the  equivalent  broad 
band  terms  of  the  susceptibilities.  Thus,  the  received  signal  can  be 
expressed  in  dB  above  1/iV/MHz. 

5.3  COMPUTER  PROGRAM  DESCRIPTION 


5.3.1  Input  Data  Required 

The  data  format  for  FTWCAP  is  designed  to  be  compatible  with  the  other- 
programs.  The  radiating  subsystem  data  and  aircraft  data  are  identical,  to 
that  used  by  ATACAP,  as  described  in  Section  3.U.I.  The  security  classifi¬ 
cation  feature  of  ATACAP  is  also  contained  in  this  program.  Further ,  th>* 
wire  bundle  data  is  identical  to  that  used  by  WTWCAP ,  as  describe i  in 
Section  U.3. 1.  The  only  additional  data  required  is  for  the  ape  rearer:  :n 
the  aircraft  skin  exposing  the  wire  bundle. 

Aperture  Data.  The  following  data  are  required  for  each  aperture 
exposing  the  given  wire  bundle: 

0  Identification  nomenclature 

0  Location  on  aircraft  (BL,  WL,  FS) 

0  Connector  locations  (BL,  WL,  FS) 

0  Bundle  lengths  to  each  aperture  (for  default  values  see  below) 

o  Aperture  length  and  width 

5.3.2  Default  Assumptions 

The  default  assumptions  in  ATACAP  and  WTWCAP  also  apply  here.  This  is 
true  for  the  aircraft,  radiating  subsystem,  and  wire  bundle  data.  An 
additional  default  assumption  applies  to  the  bundle  lengths  to  the  apertures. 
If  unknown,  the  conical  spiral  (Equation  3* LI)  Is  used  to  estimate  them, 

(Note  that  this  results  in  a  different  estimate  of  the  total  bundle  length 
for  each  aperture.) 

5.3.3  Analysis  Procedure 

The  algorithm  used  in  FTWCAP  is  summarized  in  the  flowchart  of  Figure  38. 
After  the  data  is  read,  the  default  values  calculated,  and  the  data  summer i c;; 
are  printed,  the  analysis  is  begun. 

A  transmitter-antenna  source  is  selected  from  the  data.  Using  the  same 
propagation  model  as  ATACAP,  the  RF  power  density  is  computed  at  each 
aperture  exposing  the  bundle.  The  transmission  line  equations  art*  then 
used  to  calculate  the  complex  currents  induced  in  each  susceptible  load 
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FIGURE  38  FTWCAP  ALGORITHM  FLOWCHART 
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FIGURE  38  (CONTINUED)  FTWCAP  ALGORITHM  FLOWCHART 
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FIGURE  38  (CONCLUDED)  FTWCAP  ALGORITHM  FLOWCHART 


in  the  bundle.  This  is  calculated  for  both  horizontal  and  vertical  polariza¬ 
tion,  and  the  worst-case  value  is  selected.  Included  is  any  attenuation  to 
the  RF  fields  caused  by  overall  bundle  and  individual  wire  shields.  The 
complex  signals  induced  in  the  loads  from  each  exposing  aperture  are  summed 
to  obtain  the  total  received  signal  level.  EMI  margins  for  each  load  are 
then  calculated  and  printed  when  greater  than  the  printout  limit. 

Calculations  are  made  at  the  lowest,  center,  and  highest  frequency  for 
each  transmitter  at  the  fundamental  and  at  each  significant  harmonic.  This 
procedure  is  repeated  for  each  transmitter-antenna  source  on  the  vehicle. 

5*3*4  Program  Outputs 

In  addition  to  the  input  data  summary,  FTWCAP  provides  complete  informa¬ 
tion  regarding  the  analysis.  The  information  is  divided  into  RF  power 
density  and  EMI  margin  sections.  The  propagation  data  for  each  aperture 
and  frequency  is  given  in  the  RF  power  density  section.  These  outputs 
include  the  RF  power  densities  in  watts  per  square  meter,  propagation  path 
code,  and  shading  from  wings  and  fuselage.  The  EMI  margin  section  provides 
the  information  pertinent  to  the  transmission  line  calculation.  For  each 
receptor  pin,  the  SMI  margin,  total  received  signal  (magnitude  and  phase), 
shielding  attenuation  used,  and  the  susceptibility  level  are  given.  In 
addition,  the  security  classification,  if  required,  is  printed  on  the  top 
and  bottom  of  each  page. 

5.4  APPLICATION  EXAMP  IE  AND  SAMPLE  OUTPUT 

An  actual  case  of  interference  on  an  F-4D  aircraft  provides  a  good 
example  of  an  application  of  FTWCAP.  In  this  case,  radiation  from  the  upper 
UHF  communications  antenna  was  entering  two  apertures  and  coupling  onto 
a  video  coaxial  line  of  a  receiving  system.  The  situation  is  shown  in 
Figures  39  and  40.  This  information  was  fed  into  FTWCAP  to  generate  the 
outputs  shown  in  Figures  4l  and  42.  From  this,  it  could  be  seen  that 
interference  would  occur.  (The  fix  was  to  install  an  overall  shield  on 
the  bundle.) 
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FIGURE  39  F-4D  APPLICATION  EXAMPLE 

ELECTROMAGNETIC  RELD-TO-WIRE  ANALYSIS 


FIGURE  40  SCHEMATIC  REPRESENTATION  OF  SITUATION  IN  FIGURE  39 
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FIGURE  41  SAMPLE  RF  POWER  DENSITY  FIELD-TO-WIRE  COMPUTER  OUTPUT 
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SECTION  6 


BOX-TO-BOX  COMPATIBILITY  ANALYSIS  PROGRAM  ( HTBCAP) 

6.1  BACKGROUND 

The  Increased  density  of  avionics  In  present  day  and  future  aerospace 
vehicles  has  dramatically  increased  the  problem  of  box-to-box  magnetic 
coupling.  This  applies  to  direct,  radiated  coupling.  Conducted  box-to-box 
interference  usually  involves  wire -to -wire  coupling,  which  is  covered  by 
vfiVCAP. 


An  extensive  literature  search  was  conducted,  but  very  little  previous 
work  has  been  done  in  this  area.  (Ref  8  and  18).  A  survey  of  MDC  and  USAF 
experience  indicated  that  most  direct  box-to-box  interference  problems  were 
magnetically  coupled  at  power  frequencies.  Enclosures  of  roost  avionics 
equipment  are  good  attenuators  of  RF  energy,  but  for  cost  and  weight 
considerations,  these  enclosures  are  of  low  permeability  materials  and  are 
therefore  poor  magnetic  shields. 

This  program,  therefore,  considers  only  power  frequency,  magnetic  box- 
to-box  coupling.  This  includes  the  U00  Hz  fundamental  as  well  as  the 
predominant  1200  Hz  third  harmonic  of  the  three-phase  power  systems.  The 
program  accepts  transformer  and  equipment  parameters  for  theoretical 
calculation  or  experimental  data  for  empirical  calculation.  This  applies 
both  to  source  and  susceptibility  data.  The  theoretical  and/or  empirical 
data  are  then  fitted  to  a  magnetic  dipole  mathematical  model. 

6.?  DETAILED  DESCRIPTION  OF  MODEL 

6.2.1  Source  Models 


From  physical  principles,  AC  magnetic  fields  from  the  equipment  can  be 
represented  as  magnetic  dipoles.  That  is,  the  field  source  can  be  represented 
iis  a  series  of  circulating  current  elements,  each  producing  a  magnetic  field. 
As  the  observation  distance  increases,  these  fields  merge,  and  they  appear  to 
originate  from  a  single  dipole  source.  Hence,  the  spatial  variation  of  these 
fields  can  be  modeled  using  the  magnetic  dipole  relations. 

In  most  instances,  power  transformers  in  the  equipment  are  responsible 
for  the  magnetic  fields.  The  magnetic  flux  from  the  primary  which  does  not 
reach  the  secondary  and  leaks  from  the  core  laminations  form  such  fields,  as 
Illustrated  in  Figure  Leakage  flux  accounts  for  a  large  portion  of  the 

efficiency  loss.  Although  many  other  factors  such  as  eddy  currents  and  core 
saturation  also  result  in  efficiency  losses,  for  worst-case,  all  efficiency 
loss  is  assumed  due  to  the  leakage  flux  (Ref  0  ).  The  presence  of  the  core 
material  results  in  involved  expressions  for  the  leakage  flux  near  the  trans¬ 
former.  However,  these  effects  diminish  quickly  with  distance,  and  a  dipole 
model  may  be  used. 
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FIGURE  43  LEAKAGE  FLUX  FROM  TRANSFORMER 
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Magnetic  Dipole  Model  for  Calculated  Data.  The  expression  for  the  flux 
density  from  wire  loop  carrying  current  1  along  the  loop  is  (Ref  l0) 


B  = 


M0  I  a2  N 


3/2 


where 


2  (a2  +z2) 

Q 

B  =  magnetic  flux  density  (webers/nr) 


(6.1) 


IJLq-  permeability  of  free  space 
I  =  current  in  the  loop  (amperes) 
a  =  radius  of  the  loop  (meters) 

z  =  distance  from  center  of  loop  to  observation  point  (meters) 


N  =  number  of  turns  in  loop 


As  applied  to  BTBCAP,  an  equivalent  current  is  obtained  based  on  the  trans¬ 
former  efficiency  such  that  B  represents  the  worst -case  leakage  flux. 
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It  is  given  in  terras  of  the  input  current  to  the  transformer,  Ip,  by 


I  V  1  -v 

p 


(6.2) 


where  rf  is  the  transformer  efficiency.  From  Equation  (6.2)  combined  with 
Equation  (6.1),  the  flux  density  versus  distance  from  the  transformer  can 
be  determined. 


Measured  Data.  If  measured  data  is  available,  it  can  be  fitted  to 
the  dipole  model,  and  the  model  can  be  used  to  extrapolate  the  field  at 
other  points.  Measurements  of  the  magnetic  flux  density  at  two  points  are 
required.  The  magnetic  flux  is  then  given  by 


B  (22)  =  B  (z  1 ) 


3/2 


(6.3) 


where  ae  is  the  equivalent  loop  radius, 


(6.4) 


In  the  above  zq  and  Z2  are  the  two  measurement  distances  from  the  transformer 
center,  and  B(?.^)  and  B(zo)  are  the  measured  flux  densities  in  webors/m-. 

6.2.2  Receptor  Models 

Receptors  of  low  frequency,  magnetic  interference  in  boxes  are  usually 
transformers.  Such  transformers  are  used  for  inputs  to  high  gain  audio¬ 
frequency  amplifiers,  such  as  intercom  and  servo  amplifiers.  In  the 
calculated  data  case,  a  dipole  model  similar  to  the  above  can  be  used  to 
calculate  the  interference.  If  measured  data  is  available  such  that  the 
box  is  susceptible  to  a  known  magnetic  field,  the  source  and  susceptibility 
fields  only  need  be  compared  to  obtain  the  EMI  margin. 


Calculated  Data.  The  calculated  susceptibility  is  based  on  the  frequency- 
domain  (Fourier  transformed)  version  of  Maxwell's  induction  law 

VX?=-icuB>  (6.5) 

The  induced  EMF  in  the  secondary  of  the  transformer  can  be  obtained  by 
integrating  Equation  (6.5)  over  the  loop  area.  It  is 


£=  cuBA 


(66) 
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where 


£  induced  EMF  in  loop  (volts) 
LJ  angular  frequency 
B  total  flux  density  in  the  loop 
A  -  area  of  the  loop. 


The  susceptibility  is  determined  by  assuming  the  transformer  to  be  an  M  turn 
loop  with  a  high  permeability  core.  For  worst -case,  the  medium  permeability 
is  assumed  to  be  that  of  the  core.  The  open -circuit  potential  induced  in 
the  secondary  can  be  calculated  using  Faraday's  induction  law,  and  this  is 
compared  to  the  amplifier  input  susceptibility  level  to  obtain  the  EMI  margin. 
(Note:  This  is  worst-case;  the  maximum  error  by  using  open-circuit  potential 
is  6  dB.)  This  induced  potential  is 

V,  =  co MB  (z)  pA  |6 

where  3(2)  is  the  flux  density  from  Equation  (6.1)  or  (6.3),  as  applicable, 
z.  is  the  distance  between  source  and  receptor,  jLi  is  the  receptor  core 
permeability,  and  A  is  the  receptor  core  area. 

Measured  Data.  If  data  is  available  from  MIL-STD-^l  RS01  tests,  it  can 
be  used.  In  this  test,  a  loop  of  wire  is  placed  near  the  equipment  at  a 
known  distance,  and  the  loop  current  is  increased  until  the  specification 
limit  is  obtained.  The  dipole  relation.  Equation  (6.1),  can  be  used  to 
calculate  the  equivalent  free  space  magnetic  flux  density  at  the  box  center 
location.  Thus,  the  box  is  known  to  be  susceptible  to  a  certain  equivalent 
flux  density  at  its  center.  To  obtain  the  EMI  margin  from  a  given  source. 
Equation  (6.1)  or  (6.3),  as  applicable,  is  used  to  obtain  the  flux  density 
at  the  box  center,  and  this  level  and  the  susceptibility  level  are  compared. 

6.3  COMPUTER  PROGRAM  DESCRIPTION 


6.3*1  Input  Data  Required 

The  input  data  is  easily  obtained  from  information  normally  supplied 
to  the  user.  In  some  cases,  a  number  of  options  are  available  to  aid  the 
user  in  specifying  data  in  the  form  in  which  it  is  available. 

General  Box  Data.  This  data  aids  in  identifying  the  equipment  boxes 
being  considered  in  the  run.  The  following  are  included: 

o  Subsystem  designation  or  nomenclature 

o  Humber  of  boxes  in  each  subsystem  with  respective  nomenclature 
o  Box  location  in  vehicle  (geometric  center  BL,  WL,  FS) 
o  Number  of  transformers  in  each  equipment  box 
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Individual  Transformer  Data.  This  data  specifies  parameters  associated 
with  each  transformer.  The  following  are  included: 

o  Transformer  location  with  respect  to  geometric  center  of  respective 
boxes 

o  Input  and  output  voltages  and  currents 

o  Core  material  magnetic  permeability 

o  Number  of  turns  in  primary  and  secondary  windings 

o  Susceptibility  level  of  devices  being  fed  by  receptor  transformer 
secondary 

o  Loop  area  of  coil  winding 
6.3.2  Default  Options 

The  default  options  usable  in  the  BTBCAP  are  summarized  below: 

o  Transformer  locations  are  assumed  to  be  at  the  geometric  center  of 
the  respective  boxes. 

o  When  experimental  emission  or  susceptibility  data  is  available,  the 
box  is  treated  as  a  single  transformer. 

6.3*3  Analysis  Procedure 

An  algorithm  flowchart  is  presented  in  Figure  U4.  The  input  data 
pertaining  to  each  subsystem  boxes  and  respective  transformers  are  read 
and  processed.  The  default  options  are  inserted,  and  an  input  data  summary 
is  printed. 

The  experimental  data  is  then  read  and  processed  for  emission  and 
susceptibility.  If  experimental  emission  data  is  available  then,  only  for 
emission  purposes,  the  box  is  treated  as  a  transformer.  The  same  holds  true 
for  susceptibility  data. 

The  analysis  is  begun  by  selecting  the  first  specified  receptor  and 
computing  the  SMI  coupled  to  it  from  all  source  transformers.  The  first 
source  is  selected,  and  the  methods  used  for  emission  and  receptor  data  are 
denoted.  The  EMI  margin  is  then  calculated  and  contained  in  the  output  if  it 
exceeds  the  printout  limit. 

This  process  is  continued  until  all  the  emitters  have  been  considered. 
Upon  the  completion  of  the  emitter  cycle,  the  routine  indexes  to  the  next 
receptor,  and  the  process  is  continued  for  all  receptors. 
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FIGURE  44  BTBCAP  ALGORITHM  FLOWCHART 
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FIGURE  44  (CONCLUDED)  BTBCAP  ALGORITHM  FLOWCHART 
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6.3*1+  Program  Outputs 


As  with  the  other  programs,  a  complete  summary  of  the  input  data 
including  any  default  values  is  first  printed.  Next,  the  complete  analysis 
of  each  receptor  is  printed,  specifying  the  source,  EMI  margin,  and  the 
method  of  calculation. 
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SECTION  7 


CONCLUSIONS  AND  RECOMMENDATIONS 


The  mathematical  models  and  the  computer  programs  described  in  this 
report  are  accurate,  practical,  and  useful  tools  for  establishing  intra¬ 
vehicle  electromagnetic  compatibility.  They  aid  the  user  in  determining 
and  fixing  potential  EMI  problems  on  a  vehicle  for  the  four  primary  coupling 
paths.  They  also  aid  in  establishing  the  EMC  implications  of  modifications 
to  the  vehicle  and  specification  alterations. 

The  purpose  of  these  programs  is  to  rapidly  survey  the  data  and  quickly 
indicate  potential  EMI.  As  a  result,  the  models  contain  some  approximations 
and  worst -case  assumptions.  In  a  few  instances,  the  predicted  EMI  situations 
require  further  investigation.  The  auxiliary  outputs  of  the  programs  can  be 
used  to  aid  in  these  additional  efforts.  For  example,  the  constant  wire 
inductance/capacitance  matrix  in  WTWCAP  can  be  used  in  devising  input  data 
for  a  time -domain  circuit  analysis  program  such  as  CIRCUS,  SEPTRE,  or  ECAP. 
Also,  the  ATACAP  outputs  can  be  used  to  determine  which  antenna  pairs  need 
near-field  calculations  (Ref  19) . 

This  part  presents  the  work  done  in  the  research  and  development 

phases  of  the  contract  but  not  the  verification  phase.  Consequently,  the 
models  presented  here  have  not  been  thoroughly  tested  yet  (with  the  exception 
of  ATACAP).  Where  preliminary  verification  data  was  available,  it  has  been 
presented  to  validate  the  assumptions  and  approximations.  During  the  next 
phase,  extensive  verification  tests  will  be  performed  and  the  results 
contained  in  part  II. 

The  mathematical  models  encompass  more  than  four  specific  EMI  paths. 
Because  they  are  basic  models  of  the  physical  mechanisms  of  EME  coupling, 
they  are  readily  adapted  for  other  specific  coupling  paths.  For  example,  the 
magnetic  dipole  transformer  model  from  BTBCAP  and  the  wire  circuit  model  from 
WTWCAP  could  be  combined  and  adapted  for  transformer -to -wire  EMI  calculations. 
Thus,  the  models  can  be  used  as  basic  building  blocks  for  all  combinations  of 
EMI  sources  and  receptors  including  antenna -to -box,  wire -to -box,  wire -to - 
antenna,  etc.  These  models  also  have  wide  application  to  other  areas  related 
to  intra-vehicle  EMC.  This  includes  interference  from  sources  off  of  the 
vehicle  such  as  ground-based  emitters,  other  vehicles,  EMF)  and  lightning. 
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APPENDIX 


VALIDATION  OF  LUMPED  PARAMETER  MODEL 


Consider  the  circuit  in  Figure  45. 


FIGURE  45  TRANSMISSION  LINE 

The  differential  equations  required  for  the  solution  of  this  transmission 
line  problem  are 


dx 
d  I 

—  =  —  YV  (A. 2) 

dx 

where  Z  is  the  complex  line  Impedance  per  unit  length  and  Y  is  the  complex 
shunt  admittance  per  unit  length  of  the  line.  These  equations  can  be 
solved  for  the  voltage-  across  the  load  impedance  by  differentiating 
(A. 1)  and  using  (A. 2) 


d2V 

dx2 


=  <ZY)  V 


The  solution  of  the  previous  differential  equation  is 

V  ( x )  =  AeVzY  +  Qe~*\/ZY 


(A. 3) 


(A. 4) 


Applying  the  appropriate  boundary  conditions  leads  to  an  expression 
for  the  equivalent  impedance  as  seen  from  the  generator 


ZL  cosh  (7 i)  +  Z0  sinh  I7C) 
Z0cosh  (7C)  +  ZL  sinh  (7t) 


(A. 5) 


Z0  and  7  in  the  above  are  respectively  the  characteristic 
propagation  constant  of  the  transmission  line. 


impedance 


and  the 


7  7  ZY~  (A.7) 
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Compare  this  to  the  lumped  parameter  circuit  in  Figure  46, 


z/2 
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FIGURE  46  EQUIVALENT  LUMPED  PARAMETER  CIRCUIT 

The  impedance  as  seen  from  the  generator  is 


(A. 8) 


Considering  (ZYZjJ/U  to  he  negligible  in  comparison  with  the  other  terms  in 
the  denominator  of  (A. 8) ,  this  expression  is  identical  to  (A. 5)  where  the 
hyperbolic  functions  are  expanded  to  second  order.  This  shows  that  when  the 
previous  approximations  are  valid,  the  two  formulations  are  equivalent.  For 
cases  where  the  approximation  in  (A. 8)  and  the  expansion  to  second  order  of 
(A. 5)  are  not  valid,  additional  parameters  (pi  sections)  are  added. 
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